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Many scientific discoveries and technical innovations 
were made during the tenure or NASA Grant NGR 33-010-182. 

The program involved the design, fabrication and extensive use 
of two cryogenically cooled infrared grating spectrometers on 
the NASA Lear Jet Observatory. The first spectrometer wa« 
used to me4^iure continuum sources such as dust in HIl regions, 
the Galactic Center and the thermal emission from Mars, 

Jupiter, Saturn and Venus over the 16-40ym spectral range. 

The second spectrometer had higher resolution and was used to 
measure ionic specjtral lines in 1 % regions (Sill at 18.7ym). 

It was later used extensively on NA^'A C-141 Observatory to 
make observations of ntMerous objects including HII regions, 
planetary nebulae, stars with circumstellar shells, the Galactic 
Center and extragalactic objects. It .has only recently been 
superseded by an instrument of even hi^jher resolution (built 
under NASA NGR 33-010-081) . 

In the remainder of this report the spectrometers are 
described including the major innovations and a list of the 
scientific contributions. In the final section a list of 
active participants on the project is presented. The appendix 
presents the major scientific paper published under the progi^am. 
In addition, numerous scientific papers and colloquia were 
presented. No attempt has been made to list these presentations. 

I. TECHNICAL INNOVATIONS 
A. Cryogenically Cooled Spectrometers 

Two small, liquid-helium-cGoled grating spectrometers 
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for medium resolution astronomical spectrometry in the 16-4 Oy 
region have been designed and built. They have been used in 
conjunction with NASA's 31-cm Plying Infrared Telescope on 
the Lear Jet and the 91-cm telescope of the Kuiper Airborne 
Observatory. 

A medium resolution is sufficient to allow oUservation 
of the broad spectral features associated with pressure- 
induced hydrogen absorption in the atmospheres of the Jovian 
planets. Interstellar dust found in HII regions can also be 
effectively studied with modest resolution. These instruments 
have been used to obtain spectra of the Orion Nebula and many 
other HII regions and obscured protostars. 

The design objectives for the instruments were: 

(1) a. resolution ^ ^ 30 to 200; 

(2) entrance aperture corresponding to a beam on the 
sky of 2-6-arc minutes, with the Learscope (1-3-mm entrance 
aperture) ; 

(3) cooled to liquid helium temperature to reduce the 
thermal background on the detector; 

(4) simplicity of design and operation. 

The etendue of the instrument was selected to match the 
Lear telescope (f/7.5) and the astronomical sources (<l-arc 
minute for the planets and 1-5-arc minutes for HII regions) . 

An entrance aperture somewhat larger than the image size is 
desirable because of guiding jitter caused by the telescope 
stabilization system. 

Simplicity of design and fabrication was desired 
in the instrument because of the 


near impossibility of making adjustments while it is at liquid 
helium temperatures and io maximize the reliability of the 
instrument in the severe airborne environment. Since the 
optical requirements were not severe, a single pass Ebert- 
Pastie design is used. 

B. Optical Design 

The optical layout of one of the instruments is shown in 
Figure la. The radiation enters the entrance aperture (which 
lies above the plane of the paper) and is then reflected 90° 
into the plane pf the paper by a 45° mirror. The radiation 
is collimated and then directed to the grating. Two parts of 
the diffracted beam are intercepted and sent to two separate 
detectors by concave mirrors that reduce the image size of 
the entrance aperture by a factor of 2 . 

Detector fl is an arsenic-doped silicon photoconductor 
which is normally used in second order for wavelengths of 16- 
23ii. Detector #2 waii a gallium-doped germanium photoconductor 
(it was later replaced by a zinc-doped germanium photoconductor) 
which is used in first order for wavelengths of 20-40vi. Thus 
for each grating position a short wavelength and a long wave- 
length are observed simultaneously; this procedure shortens 
the time necessary to measure a given spectrum and allows one 
to determine the extent of guiding jitter. 

Detector #1 can be used in first, order, if desired. Then, 
it will respond to radiation from wavelengths of 16-24iJ. 

The spectral orders are separated by an interference 
filter for the short wavelength detector (#1) and by the two 
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CaS ’2 restrahlep, mirrors for the long wavelength detector. 

The only • transmission element in the system is the KRS-5 

Dewar window. 

The grating is a sitandard Bausch and Lomb replica; it 
has 21 grooves per millimeter and is blazed for 28.56)j, The 
ruled area is 26x26 mm, which is not quite fully illuminated. 
The focal length of the primary mirror is 125 mm, which gives 
a dispersion of 0.38)j/mm in first order and a nominal resolu- 
tion of 1. Op in first order and 0.5p in second order with our 
2.7-mm diameter entrance aperture. 

The optical layout of the second instrument is shown in 
Figure lb. The entrance aperture and collimating optics are 
I similar to those described above. The dispersed beam is 

filtered by restrahlen reflection plate and an interference 

% 

filter. The instrument was firs.;; flown with a four-element 
array of cooper-doped germanium detectors. Although the in- 
strment performed well, the detectors had 2 or 3 times worse 
NEP's than corresponding Si; As elements. The performance was 
dramatically improved by the installation of a ten-element 
array of doped-silicon detectors (Si;Sb). 

C. * Mechanical ■Considerations 

The spectrometer itself is 125 mm (5 in.) in diameter 
and 59.5 mm (2-11/32”) tall and is mounted on the OFHC copper 
baseplace of a liquid helium can inside a helium Dewar as 
shown in Figure 2. The spectrometer is constructed entirely 
* of alumin^Im (6061-T5) except for the detector cavities Which 

are of OFHC copper. 




The Ge:Ga (GesZn) detector is heat sunk directly to the copper 
baseplate of the Dewar,- while,; the silicon detectors are heat 
sunk to the baseplace of the spectrometer. 

The primary mirror and the grating are mounted in spring 
loaded holders with copper cold straps to the bcdf'y of the 
spectrometer . 

All the small mirrors are attached to their mounting 
posts with either rubber cement or epoxy. These mirrors are 
supported only at one point to minimize distortion of their 
surfaces due to the differential thermal expansion between 
aluminum and glass or CaF 2 . In many years of experience with 
more than IQO cool«down cycles, none of these joints has 
failed. No problems have been encountered with any of the 
optical components which can be associated with the helium 
temperature Operation . 

The grating carrier is attached to shafts which ride in 
ball bearings at the top and bottom of the spectrometer. On 
the bottom is a sector of a worm gear mounted on the grating 
shaft. A worm is also mounted on the bottom plate and meshes 
with the gear. The worm is connected via two flexible bellows 
couplings and a fiberglass shaft to the outside of the Dewar 
through a vacuum tight 0-ring fitting. Once outside, the 
grating driveshaft is connected to a Superior Electric step- 
ping motor by 1:1 bevel gears. A shaft encoder and a rotation 
counter on the worm shaft are used to monitor the grating 
position. / 


Spect ?' ometer Performance - Leau: Jet 30 -cm Telescope 

• 10 -detector array spectrometer 

XX16-30 AX » 0.5U Si:Sb detectors 

Beam size « 2.7' 

AH « 3.6xl0“‘*cm^ ster 

In September-October, 1977, we observed VY CMa on several 
flights. From our flux calibration we estimate P 23 (VY CMa) 
tt5 X 10"l5w/cm2y which corresponds to F^(23y, VY CMa) « 8800Jy. 
From an extensive series of noise measurements during several 
flights, we estimate an rms S/N ratio of approximately 22 for 
one second of integration time at wavelengths between atmo- 
spheric water vapor absorption. 

If we define; 

N®‘system ^ . (1) 

and 

NES’Dgystem * s7W 

where t = integration time, P - power from the source in the 
detector bandpass incident on the telescope, F^ is the flux 
density of the source, S/N is the estimated rms signal-to- 
noise ratio. Then the above data imply a system performance 
at 23 ]i wavelength of: 

system = 8 x lO-'V^ 

and 

*’^^°system ” 

for each 0.5y bandpass. The individual dsstectors are matched 
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to within ±10% of this perfoimance, 

Poir the observed in-flight detector performance, we esti- 
mate the spectrometer transmission to be ^20%, the beam-split- 
ter reflectivity to be a90%, and the chopper efficiency to be 
«90%. Then the observed detector performance in flight is 
approximately 


“ 1-3X10-‘V/W 

for each of the 10 detectors. 

•• Two-Detector Spectrometer 

a 


XX16-23 AX = 0.5p Si:As detector 

XX21-38 AX l.Ou Ge:Zn detector 

beam size » 3.3' 

AJ] « 5.5 X 10"*^ cm^ ster 

From measurements made on VY CMa in November 1976 , we 
estimate, using equations (1) and (2) , and the techniques out- 
lined above, in flight performances of 


Si:As @ 19.3 

'“'’system “ “ 

"''''“system “ 

"'“detector “ 


Zn;Ge Q 31.2 
8 X 10"! %//Hz 
3900 Jy//Hz 

1.3 X 10-13^//H^ 


E. Servo-Controlled Chopper Systems 

The early choppers used on the Lear Jet employed solenoids 
and mechanical limit stops. As a result no adjustments could 
be made of the chopper during flight to either change the 
chopper throw or its offset. Because of the mechanical stops, 
the chopper was very microphonic and required frequent adjust- 
ment. These early choppers were replaced by a servo-controlled 
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chopper in 1975. The chopper was designed at -Cornell and 
fabricated in the An,ies model shop. The electronic drive unit 
was designed and built at Cornell. The system was optimized 
at Cornell and field-tested on the Lear Jet. It was subse- 
quently turned over to the Ames personnel and has been maintained 
by them since 1976., Since then it has been the standard 
chopper for all Learscope observers who use choppers. 

lar. SCIENTIFIC RESULTS 

Numerous scientific results have been achieved using these 
spectrometers on the Lear Jet. In this section a very brief 
list of these is presented. The appendix includes a much 
more complete summary of our results. 

Planetary Objects: 

Jupiter - Determination at the He/H abundance ratio. 

Saturn - Same aiS for Jupiter 

Mars - Determination of continuum spectrum. 

Venus - Detection of spectral features due to H 2 SO 4 in 
the atmosphere. 

Galactic HII Regions ; 

1) Determination of the shape of the ”20ym silicate 
feature” in a large number of sources. 

2) Measurement and interpretation of the Sill 18.7ym 
spectral line in nearly a dozen HII regions. 

Stars : Determination of the long wavelength distribution 

of emission from VY Canis Majoris. 

Galactic Center: Measurements of the 16 to 30ym continuum 

emission from the Galactic Center. 
Determination of the dust extinction 
toward the Galactic Center. 
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III. PERSONNEL 

The success of this program is dp to the efforts of many 
people. At Cornell these included three students who received 
their Ph.D.'s for research conducted on this program. Dr. W. J. 
Forrest was a post-doctoral research associate during much of 
the program and was the central figure in much of the work. 
Extraordinary technical support was provided by G. Stasavage, 

G. Gull and D. Briotta. The efforts of the Ames personnel were 
crucial to the success of the program. We particularly appre- 
ciate the assistance of the Lear Jet Manager/ Robert Mason, 
and his staff. We also appreciate the Lear Jet pilots who 
spent many uncomfortable and often very cold hours in the middle 
of many nights helping us "get the data." 

PH.D. THESES 

David Francis Schaak, 1975, "Infrared Astronomical Spectroscopy 
from High Altitude Aircraft." 

Robert Alan Reed, 1976, "The 16-to 40-Micron Spectra of Jupiter 
and Venus . " 

John Francis McCarthy, 1980, "Airborne Infrared Spectroscopy of 
Ionized Hydrogen Regions and the Galactic Center. " 
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20 TO 40 MICRON SPECTROSCOPY OF THE ORION NEBULA 

J. R. HOX7CK,* D. F. SCHAACK, AND R. A. Reed 
Center for Rediophyslcs and Space Research, Cornell University 
Rtctivtd 1974 JiUy 22 

ABSTRACT 

We present moderete-resoluticn (1.0 ft) spectra of the Orion Nebula. These observations are the first spectroscopic 
observations of the object between 20 and 40 ft. A comparison of our data with that of Harper, Solfer and Hudson, 
and Harvey tt al, is presented. We find that the best fit to all of the data is provided by a model in which the dust 
within the molecular cloud has a clumpy distribution. The individual clumps have effective grain temperatures of 
100* K and are optically thick for wavelengths shorter than 400 ft, The total dust moss required by this model is less 
than 20 Mq, 

Subjici htaJints: infrared — intentellar matter — Orion Nebula 


I. mTEODUCTtON 

Two types of extended infrared sources have been 
detected in the Orion Nebula. The fint type shows a 
strong 10-M emission feature thought to result from 
circumstellar shells of silicate dust. The Ney-Allen 
(1969) object is of this type. The second class of objects 
shows the same feature in absorption (e.g,, the Klein* 
mann-Low Nebula). These objects are thought to be 
‘‘stars” embedded in much larger dust clouds. The 
former objects seem to cluster in the direction of the 
Trapezium while the latter are centered about the posi- 
tion of the molecular cloud (Zuckerman 1973). 

As first shown by Harvey el d. (1974), the bulk of 
the submillimeter radiation also comes from the direc- 
tion of the molecular cloud. At 1 mm the apparent size 
of the cloud is ^^2'. We present new observational data 
between 20 and 40 fi. Several models are presented in 
an attempt to account for all of the obser/ations from 
20 ju to 1 mm. 

Analysis of the 28.2-)u line of molecular hydrogen 
will be presented in another paper. 

rr. OBSERVATIONS 

The Orion nebula was observed on the nights of 1973 
November 14-21 using the 31-cm telescope mounted 
on the NASA Lear jet. This telescope is very similar to 
the one developed by Low. The aircraft was flown to 
an altitude of 14 km (45,000 feet). At this altitude about 
3 to 10 prccipitable microns of water per air mass re- 
main above the observer. The telescope viewed the 
source at elevation angles between 14° and 28° (4.1 
to 2.1 air mass). The spectral scans were made using a 
iiquid-helium-cooled grating spectrometer (Schaack 
and Houck 1974). This instrument employs a 12.5-cm 
focal length Ebert-Fastie spectrometer with two de- 
tectors in the focal plane. A Ge:Cu photoconductor 
scans the 16-28 m hand while a Ge: Ga photoconductor 
covers 20-40 fi. The resolutions of the two channels are 
0.5 and 1.0 ft, respectively. We report here only the 
results of the 20-40 ft channel. The instrument operates 
at f/7.5 and has a 2.7-nun entrance aperture corre- 

* Alfred P. Sloan Foundation Fellow. 


spending to 4(7 on the sky. A wobbling secondary 
cnopper was used. 'Ihe beam separation was 15', The 
chopping frequency was 48 Hz. 

The instrumental response was determined by nor- 
malizing spectra of the Moon and Mars as if they were 
blackbody radiators at 350° and 220° K, r«*spectively. 
The instnunental profile determined in this way was 
used to noimalize the Orion data. The resulting spec- 
trum is shown in figure 1. The spectrum was normalized 
to the 30-45 fx broad-band photometric point deter- 
mined by Harper (1974). This is consistent with an 
absolute calibration of the system by intercomparing 
Lear-jet and ground-based observations of Jupiter. It is 
also consistent with the low-spatial-resolution 20-/* 
ground-based map of Orion obtained by Lemke, Low, 
and Thum (1974). The major noise source was tracking 
error associated with the pointing system of the Lear- 
scope. 

m. DISCUSSION 

I 

Figure 2 shows representative far-infrared obsewa- 
tions of Orion. The dashed curve is the result of a model 
in which the source is a uniform 2' diameter cloud as 
suggested by the observations of Harvey ei d. It has 
been assumed that the grain emissivity is proportional 
to the reciprocal wavelength squared. This is the form 



Wavaiength XC/il 


Fio. 1. — Spectium for the Orion region from 20 to 40 ft. The 
beam diameter was 417. The curves are of blackbody radiation at 
various temperatures. 
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which results from using the optical constants for lunar 
material and the scattering theory for small particles 
(^acke and Thompson 1973), Hoyle and Wickrama- 
singhe (1962) have also suggested this frequency de- 
pendence for graphite grains. The model was fitted to 
the data joints at 30 and 1000 fi, This two-parameter 
fit determined the grain temperature T, and the wave- 
length at which the cloud becomes optically thick, Xr; 

Fobs = ■5(X, rB)CNEB[l — e-xp — (Xr/X)’] , (1) 

where Fobs is the observed flux, B(\, T) the Planck 
function, Oneb the solid angle of the nebula (Uneb = 
2,1 X 10“’’ sr for a 2' disk). For this model Tf dS” K 
and Xr = 180 p.. 

The solid curve results from a model in which the 
emissivity again is proportional to In this case, 
however, an attempt was made to fit the shape of the 
observed data disregarding temporarily the absolute 
flux comparison. Under these conditions the best fit is 
obtained by Xr => 400 /u with Tf =* 100® K. A direct 
application of equation (1) predicts 10 times more flux 
at each wavelength than is observed. Strong clumping 
within the source would of course reduce the flux below 
that predicted by equation (1), In this case one would 
conclude that only 10 percent of the 2' beam was filled 
by sources each with Te = 100® K and Xr ™ 400 p. 


Rjekc, Low, and Kleininann (1973) have mapped the 
region around the Kleinmann-Low nebula at 5, 10, and 
20 p with very high spatial resolution (5"-10"). They 
find that the region Is divided into at least five separate 
sources, These clearly separated sources arc too faint to 
account for the flux Observed between 30 and 45 p by 
Harper (1974) from the Lear jet by a factor of 10 to 20, 
assuming the small sources have the same spectral 
shape as the entire region. Riekc et al, also find evidence 
at 10 p for an underlying smooth (or unresolved distri- 
bution of sources) flux distribution. Even including this 
radiation, the ground-based observations detect only 
10-20 percent of the large beam flux observed from air- 
craft. This discrepancy is probably due to the self- 
cancellation inherent in beam-switched observation of 
extended sources. Typical beam throws for airborne 
observations arc 10', whereas they are is seldom larger 
than 1' for ground-based observations. Therefore, 
ground-based observations support a clumpy model, 
but we are unable to make a clear connection between 
ground-based and airborne observations. 

The “uniform cloud model" is pressed to explain any 
small-scale features, If the cloud has Xy « 180 p, then 
the optical depths at 20, 10 and 5 p are 15, 20, and 15, 
respectively, if the data of Knacke et of. are used. 
Under these conditions the sources seen by Rieke et aU 
would have to be very near the surface of the cloud to 
be seen at all. The sources would be confined to ^5 per- 
cent of the volume of the cloud. 

A third model involves a “uniform cloud” with 10 
percent of the area of the 2' source seen by Harvey et al, 
with Tf - 100® K and Xr 400 p. This model would 
result in an enormously bright object 40^ in diameter 
that is inconsistent with ground-based surveys* 

A fourth model is one in which Q varies as X“® for 
X > 400 p and Q » constant for X < 400 p. Such a 
behavior would be quite unlike that of lunar material. 
However, the results of Knacke et al, are derived from 
room temperature measurements of materials that have 
undergone sufficient heating to acquire recognizable 
mineralogical forms. For the Orion cloud the tempera- 
ture is certainly much less than room temperature. 
Further, it is quite likely that the grains have much less 
regular structure than lunar fines. Hoyle and Wick- 
ramasinghe (1969) suggested that the impurities (in- 
cluding structural imperfections) would lead to en- 
hanced emission at long wavelengths. One would there- 
fore invoke enhanced emission due to a variety of 
impurities to maintain^ constant with r as 0.1 from 
20 to 400 p. This degree of enhancement would be 
possible based on the detailed analysis of impurity 
mode emission bj’ Field (1969). At wavelengths longer 
than 400 p one would still expect a X“® dependence as 
shown by the arguments of Gezari, Joyce, and Simon 
(1973) and Gezari et al, (1974), 

The mass of dust involved is best detennined by 
fitting the parameters of the source and dust at 1 mm 
where the optical depth is small for all models; 
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where R is the distance to the nebula, 460 pcj Q is the 
ratio of effective emitting cross-section to geometrical 
cross-section; and a is the grain radius (taken to be 
0.2 fi). The value of Q/a ranges between 0.3 cm”* 
(Knacke and Thomson 1973) to 400 cm"* (Werner and 
Salpeter 1969). The higher value corresponds to the 
ma.'dmum allowable value based on the inclusion of 
impurity oscillators. 

For the clumpy model the mass of the dust is 0.02- 
20 Me depending on the choice of Q/a. The uniform 


disk model requires 50 percent more dust because of the 
lower temperature. Using dynamical considerations, 
Zuckerman (1973) has estimated the total molecular 
hydrogen mass to be < 10* Mo, The above estimates are 
tlierefore consistent with cosmic abund?/nces. 

We wish to thank J. L. Pipher and B. T. Soifer for 
many helpful discussions. We are also grateful for the 
efforts of the NASA Ames Airborne Science Office. This 
work was supported by NASA grant NCR 33-010-182. 
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Jupiter; Its Infrared Spectrum from 16 to 40 Micrometers 

Abstract. Spectral measurements of the thermal radiation from Jupiter In the band 
from 16 to 40 micrometers are analyzed under the assumption that pressure-broadened 
molecular hydrogen transitions are responsible for the bulk of the Infrared opacity over 
most of this spectral Interval. Both the vertical pressure-temperature profile and the mo- 
lecular hydrogen mixing ratio are determined, The derived value*of the molecular hydro- 
gen mixing ratio, 0.89 ±0,1 1, Is consistent with the solar value of 0.86. 


The abundance of H, and He in the Jo- 
vian atmosphere has a direct influence on a 
number or astronomical problems. Jupi- 
ter's low atmospheric temperature and 
great mass prevent even the lightest atoms 
from escaping from the top of the atmo- 
s|.here. Therefore, Jupiter as a whole is a 
sample of the elemental abundance at the 
time of the formation of the planet. Mod- 
em theories of explosive nucj ^synthesis 
show that most of the He pre»..itdy in the 
universe was formed during the Big Bang 
(/). Further, it has been shown that the rel- 
ative abundance of H{ and He depends 
strongly on the temperature and density 
during the early stages of the evolution. To 
the extent that the Jovian atmosphere is 
representative of the plan<?( as a whole, a 
measure of the H, mixing ratio oh, - 
.\(H})/(,V(H,) + ^{He)|, where N{x) is 
the number density of ,r in the atmosphere, 
is useful for the determination of the con- 
ditions during the early stages of the Big 
Bang. There are, however, several effects 
which may systematically distort the at- 
mospheric value of OH,. 

Opik has suggested a model for the for- 
mation of Jupiter in which first “hydrogen 
snow" collects to form a core and then a 
He-rich atmosphere is captured (2). Salpe- 
ter has suggested that there may be inter- 
nal differentiation with the He sinking to- 
ward the center of the planet, resulting in a 
Hi-rich atmosphere {J). 

We observed Jupiter on the nights of 14 
November and 16 November 1973 and 21 
January 1974, using a 31-cm telescope 
mount^ on the National Aeronautics and 
Space Administration Ames Research 
Center Lear jet. The aircraft was flown to 
an altitude of 14 km (45,000 feet). At this 
altitude about 3 to 10 precipitable microm- 
eters of water remain above the observer. 
The telescope viewed the source at eleva- 
tion angles between 14* and 28*. The spec- 
tral scans were made with two different 
grating spectrometers cooled to liquid 
helium temperatures {4), These instru- 
ments employ an Ebert-Fastie spectrom- 
eter (12.5-cm focal length) with two de- 
tectors in the focal plane. Instrument I has 
a Ge:Cu photoconductor to scan the 16- 
to 28- am band, and a Ge : Ga photocon- 
ductor covers the range from 20 to 40 Min. 


The resolution of the two channels is O.S 
and 1,0 *Mm, respectively. Instrument 2 has 
two Ge : Cu photoconductors and covers 
the range from 16 to 28 Mm with a 2.7-mm 
entrance aperture corresponding to 4.7 
minutes of arc, Jupiter had a diameter of 
about 0,5 minute of arc during the ob- 
serving periods. The chopping frequency 
wds 48 hertz. 

We determined the instrumental re- 
sponse by normalizing spectra of the moon 
and Mars as if they were blackbody radi- 
ators at 350* and 240*K, respectively, The 
instrumental profile determined in this way 
is consistent with the laboratory profile 
and the absorption expected by the atmo- 
sphere. It was used to normalize the Jupi- 
ter data. The resulting spectrum (Fig. la) 
shows the presence of three absorption fea- 
tures at 18, 23.5, and 28 Mm. The first and 
last are. respectively, the y ■ I and y - 0 
rotational transitions of H.. The 23.5- Mm 
feature shown in Fig. la may be due to 
sulfur, silicate dust, or complex hydro- 
carbons in the atmosphere (i). However, 
recent data of greater resolution and sig- 
nai’ to-noise ratio indicates that the bright- 
ness temperature falls smoothly from a 
maximum at 21 Mm to a broad minimum 
around 28 Mm. The regions of the spec- 
trum used to determine the ratio of Hj to 
He are substantially unchanged. However, 
the reality of the 23.5 Mm is somewhat in 
doubt. A discussion of ail the observational 
data will be given by Pollack et al, (6), 

The observed spectrum of Jupiter con- 
tains information about both a portion of 
its vertical temperature structure and its 
ratio of He to H,. Before describing our 
numerical method for deriving this infor- 
mation, we describe the physical connec- 
tion between these quantities and the ob- 
served spectrum. The brightness tempera- 
ture found at a given wavelength is approx- 
imately equal to the value of the physical 
temperature in the Jovian atmosphere at 
an optical depth of unity. As the wave- 
length changes, the altitude at which the 
optical depth is unity also varies. Thus, 
spectral observations over a range of wave- 
Jengths provide information on the temper- 
ature conditions over a corresponding 
range of altitudes in the atmosphere. For 
the spectral band measured. Hj is the prin- 


cipal source of opacity. Jovian at- 
mosphere is prolKd from pressures of 
about 0.15 t'} 0.6 atm, a region which in- 
cludes a temperature minimum and the top 
ofthe convection zone. 

In carrying out the analysis described 
below, we excluded data close to the 23,5- 
Mm feature. In addition to opacity due to 
the rotational and translational transitions 
of Hj. we allowed for opacity due to NHj, 
which is important only at the long- 
wavelength edge of our data. We carried 
out calculations both without and with ai) 
optically thick NH, cloud present in the 
zone regions of Jupiter, located at the level 
of the atmosphere where NH, saturation 
occurs. The NH, clouds were assumed to 
behave as blackbody emitters. 

In accord with radio temperature mea- 
surements we assumed that the temper- 
ature gradient equaled the adiabatic value 
at pressures greater than the bottom alti- 
tude boundary of the region we sense (7), 
Our derived temperature profiles are con- 
sistent with this assumption. Above our 
top altitude boundary we made use of the 
temperature gradients derived by Orton (8) 
from an analysis of measurements in the 
region of the 7.7- mhi C H , band. 

Our ability to gain information about 
the ratio of He to H, follows from the de- 
pendence of the pressure-induced transi- 
tions of H, on litis ratio. This dependence 
is illustrated in Fig. lb, which shows the 
values of the absorption coefficient of H, 
when the rotational transitions occur in the 



Fig. I . (a) The brighincss temperature of Jupiter 
for wavelengths from 16 to 40 Mffl. (b) The pres- 
sure-indiKed absorption coeflicienu for H.-He 
collisions. The curves shown represent the sums 
of both translational and routional transitions. 
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Fig. 2, A plot or the derived pressure-tempef 
ature profile for the Jovian atmosphere. The 
pressure level most directly measured by various 
observational techniques is shown. 


presence of a nearby molecule of H, (H,- 
H j), and in the presence of a nearby He 
atom <H,-He). We see from Fig. lb that 
the absorption coefficient shows a different 
spectral variation for the (Hj*He) case 
than for the case. 

The numerical method used to obtain 
the desired information from the observed 
spectrum of Jupiter begins with an as* 
sumed trial value of oh,. For this choice of 
the mixing ratio, the observed spectrum is 
inverted to determine the pressure*temper- 
ature structure by means of an iteration 
technique developed by Smith (9) and ap> 
plied, to planetary atmospheres by Ohring 
(10). This method is similar to the method 
developed by Encrenaz and Gautier (11) 
and Gautier and Grossmann (12), It in- 
volves starting with a trial temperature 
profile, computing a predicted spectrum 
with this profile, and correcting the profile 
on the basis of the difference between the 
observed and calculated fiux values, the 
fiux residuals. Tlie improved guess is used 
to repeat this process and, after a small 
number of iterations, the root-mean- 
square of the fltix residuals approaches an 
asymptotic value. This procedure is re- 
peated for other trial values of oh,. An es- 
timate of OH, was obtained from the loca- 
tions of the minimum in the value of the 
flux residual as a function of mixing ratio. 

Figure 2 shows the vertical temperature 
structure of Jupiter determined from our 


spectrum for an, values of 0,8 and 1 , 0 , 
These results refer to models containing no 
NH, cloud opacity. However, very similar 
curves were obtained for measurements 
with optically thick clouds. Our results 
have been smoothly joined with the tem- 
perature gradients* found above and beiow 
our region of sensing. Figure 2 also in- 
dicates the location of the regions probed 
in our measurement, at radio wavelengths, 
and within the 7,7- Mm CH 4 band. The tem- 
perature values found from these various 
determinations are quite consistent near 
their boundaries. The temperature gradi- 
ent becomes adiabatic at the higher-pfes- 
sure domain of our sensing region, which is 
consistent with the temperature lapse rate 
implied by the radio results, We find that 
the lapse rate first reaches the adiabatic 
value and therefore that the convection 
zone begins at approximately 0,4 atm. This 
result is in good agreement with the theo- 
retical predictions of Pollack and Ohring 
(13), which were based on a radiative equi- 
librium model. Our measurements also 
suggest the presence of a temperature min- 
imum and the start of an inversion layer at 
the lower-pressure levels of our region. 
This aspect of the temperature profile is 
suggested by the small displacement of the 
observed minimum in our brightness tem- 
perature spectrum from the center of the 
/ ■ I H] rotational transitions. These re- 
sults are in accord with Orton’s analysis 
( 8 ) of measurements made in a spectral re- 
gion containing the 7.7- nm fundamental 
of CH 4 . The value of our temperature min- 
imum (II3*K) is also in good agreement 
with' Orton’s value of about I I 8 *K. 

The temperature profiles shown in Fig. 2 
are in conflict with the structure found 
from a preliminary analysis of the Pioneer 
10 S-banti' occultation experiment (14), 
The results from the occultation experi- 
ment indicate a much warmer atmosphere 
than our results at pressures above several 
millibars. Our infrared observations could 
be reconciled with the S-band data if we 
postulate the presence of an optically thick 
cloud near the I -mbar region. However, in 
this case the pressure-induced opacity of 
H] would be negligible and we would not 
expect to detect the J ■ 0 and 7 ■ I ro- 
tational transitions in our spectrum. The 
presence of these features in our spectrum 
suggests that further study of the S-band 
occultation results is needed. These con- 
clusions are also supported by the con- 
sistency of our temperature profile with 


values obtained at shorter infrared 
wavelengths and in the radio domain, as 
discussed above. ^ 

For models containing no NH, cloud 
opacity in the zones of Jupiter, we find that 
the OH, is 0.85 ^ 0,07, whereas we obtain a 
Value of 1.0 t'lat for the cases Involving an 
optically thick NH, cloud. These error 
bars reflect only the formal random errors 
of our results. The cloud-lrae models had a 
smaller value for the minimum fractional 
flux residual than the set Inco'hiprating a 
thick cloud (0.039 versus 0,043), Com-’ 
bining the above two estimates, we con- 
clude that n-;;, equals 0,89 ;S: 0.1 1, This 
range of values encompasses the- solar 
value of 0 . 86 . 
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ABSTRACT 

Observations of the 16-40 m spectra of the Trapezium region and the Kleinmann-Low nebula in 
the Orion Nebula with a resolution VAX — 30-40 are reported. The Trapezium spectrum shows a 
single broad emission feature peaking at 18.3 this is interpreted as the O-Si-0 bending^mode 
resonance in silicate grains. The spectrum is compared with a simple model of small heated silicate 
grains. The Kleinmann-Low Nebula shows a continuum with the color temperature approximately 
equal to brightness tei^erature {'--83 K). The weakness of the 20 m absorption feature is attributed 
to radiation transfer effects in the molecular cloud. .A decrease in brightness temperature to longer 
wavelengths indicates that the center of the molecular cloud is being observed at 28 n. The cloud 
mass and visual optical depth derived using this information are in good agreement with molecular 
determination while the derived central density is an order of magnitude larger. 

Subject headings: infrared: spectra — nebulae; Orion Nebula 


I. JJITSODUCTION 

.A general description of the Orion Nebula sources 
and the importance of 20 u spectroscopy is given in the 
accompanying Letter (Forrest and Soifer 1976, herein 
referred to as FS). The present spectra are important 
in e.xtending the wavelength coverage beyond 24 u- 
With the entire 16-40 u spectrum a more complete 
picture of the emissivity’ properties and temperature 
and optical depth of the dust in these sources is formed. 
In addition, the present spectra should aid in the 
interpretation of the longer wavelength observations of 
these sources. 

II. INSTRUJIENTATION 

The observations were made with a helium-cooled 
Ebert-Fastie spectrometer (Houck et at. 1976). An 
As: Si detector was used for the 16-23 fi band. The 
20-40 u band used a Ge:Ga photoconductor built bv 
W. J. Moore. The short- and long-wavelength band- 
passes were 0.3 jx and 1.2 n FWHM, respectively. The 
spectrometer was used in conjunction with the *91 cm 
telescope of the NAS.A Gerard P. Kuiper Airborne 
Observatory. The beam size was 23", and the beam 
throw was '—90" in an E-W direction. The C141 
aircraft was flown at an altitude of 12.7 km (41,000 
feet). The total observing time available for the 
observations reported here was hours (plus 1 hour's 
observation of Mars). 

in. OBSERVATIONS 

The Kleinmann-Low Nebula (Kleinmann and Low 
1967, herein referred to as KL) was observed on two 
flights, 1973 October 29-30 and October 30-31, The 
Trapezium region (Ney and .Allen 1969) was observed 
only on the second flight. Mars was used as a calibration 


source and was observed on both flights. An on-board 
water-vapor column-density mete», indicated 7-15 n 
of precipitable water vapor above the aircraft during 
the first night, and the Martian spectmm on the 
second night was in good agreement with the first 
flight spectrum, with no indication of differing terres- 
trial water vapor, The typical night to night repro- 
ducibility of the KL spectra was '-*'5 percent. The 
observed signals were converted to flux levels by using 
Mars as a calibration source. Its ephemeris diameter 
and a blackbody temperature of 240 K were assumed,' 
.A small «-20 percent wavelength-dependent correction 
for \ u was used to account for COi absorption 
in the Martian atmosphere. This correction was de- 
termined by comparing spectra of Mars and the Moon 
using the spectrometer on NASA’s 31 cm Lear Jet 
telescope. Figure 1 shows the derived spectra under 
these assumptions. 

.Also included on Figure 1 are ground-based broad- 
band measurements of these sources made by other 
observers with comparable beam sizes. The agreement 
with the 18 and 20 m measurements of Ney, Strecker, 
and Gehrz (1973, 25" beam) and Ney and Allen (1969, 
26" beam) is fair. The early 22 ji measurement of KL 
by Kleinmann and Low (1967) seems too_ low — since 
these observers used a 30" beam, the disagreement 
cannot be ascribed to beam -size effects. The 34_ fi 
measurements by Low, Rieke, and Armstrong (1973) 
with a 23" beam also seem systematically low compared 
to the present observations. Part of this discrepancy 

' Wright (1976) has proposed a model for the thermal emission 
from the whole disk or Mars. .At the time of these observations, 
the model gives a Martian brightness temperature of 240 K at 
—30 ft but a slightly higher (—260 K) color temperature from 
20 to 40 ft. Th& would not materially aflect the conclusions 
reached here. 
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Fio. I,— The 16-38 spectra of the Klelnmann-Low (KL) nebula and the Trapeefum region In the Orion Nebula. The smooth curve 
represents an 8S K blackbody filling the 25* beam, For the KL spectrum the statistical errors may be estimated from the scatter of adja- 
cent points. For the Trapezium spectrum the sutistical errors from 16 to 23 (ehamsel 1) can be estimated from the agreement of the 
second spectrum obtained (r/dsrd circles) with the first (open circles); for channel 2 (20-38 sj. open souares), data over approximately 2 ft 
passbands have been averaged and error bars representing 1 <r of the mean are shown. Broad-band data from Ney and Allen ( 1969, open 
diamonds), Ney el al, (1973, open Iriangles), Low el al, filled diamonds), and Kleinmann and Low (1967, solid Iriangle) are included. 


for the Trapezium could be due to the smaller beam 
throw ('"-'45") these authors employed. The discrepan- 
cies at 34 M are believed to be larger than the systematic 
errors in the present e.'tperiment. 

IV. TITE TSAPEZmi 

The spectrum of the Trapezium (Fig. 1) shows a 
broad emission feature peaking at ~18.5 n with flux 
decreasing out to the longest wavelengths. No emission 
or absorption lines are evident to within '"-10 percent 
of the continuum level. This spectrum agrees quite 
well in shape with the spectrum presented in the 
accompanying paper (FS)j the only real difference is 
that the present spectrum appears to decrease slightly 
less rapidly longw’ard of 19 n. The emission feature is 
narrower than a blackbody; therefore, it indicates a 
resonance in the emissivity of the dust grains near 
18.5 ft. This is in qualitative agreement witli the 
expected behavior of silicate grains. 

Figure 2 shows a plot of both the long (16-40 m) and 
short (8-13 /i) wavelength spectra of the Trapezium, 
region. The short wavelength data represent the spec- 
trum reported by Forrest, Gillett, and Stein (1975) 
scaled to a 25" beam; the adopted 10 ft flux level agrees 
with the measurements of Gehrz, Hackwell, and Smith 
(1975) with a 26" beam. The dashed curve is a model of 
160 K silicate grains fitted to the data. The emissivity 
as a function of wavelength is for Moon rock No. 14321 
(Perry el al. 1972) in the form of 0.2 n grains taken 



Flo. 2 . — .\ simple silicate model (dashed curve) fit to the Trape- 
zium data ai described in the text. The open squares and closed 
circles are the present data from Fig. 1, and the open circles are 
the present dau with 3% of the flux from KL (Fig. 1) subtracted. 
The solid line represents the 10 u spectrum repotted by Forrest 
et ai. (1975), 
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from Knacke and Thomson (1973). The optical depth 
of this model is ~-0.3 X 10"* at 20 n, 

The actual temperature of the dust radiating at 20 n 
is difficult to deduce, as it is related to the relative 
strengths of the 10 ii and 20 m emission features. In 
addition, Dav (1976) has shown that the relative 
strengths of these bands can vary with grain tempera* 
tures, the 20 m band becoming stronger at lower 
temperatures, For ratios irom i to I, the dust 

temperature would range from ISO K to 200 iv. These 
temperatures are significantly lower than the 250 K 
temperature deduced from the 10 m emission profile 
(Gillett W al. 1975), This discrepancy could possibly 
be e.xplained as an effect of a temperature distribution 
within the Trapezium, or it may indicate that the cosmic 
silicates have a larger 20-10 ju emissivity ratio than 
has been reported from laboratory measurements. Since 
Ney, Strecher, and Gehrz (1973) found only a small 
(~2S K) range in the 10-20 m color temperature for 
several regions in the Trapezium, the latter effect would 
seem to be more important, 

It is seen from Figure 2 that the simple model gives 
a fair representation of the overall shape of the observed 
spectrum from 8 to 20 m. Longward of 20 there Is 
increasing disagreement in that the observed spectrum 
falls less rapidly than the single*temperature silicate 
model. There are Several possible explanations for this 
discrepancy. The grain emissivity could be falling less 
rapidly than the lunar silicates for which t>, ~ X"’ 
here; for 160 K grains it requires an emissivity approxi- 
mately constant from 25 to 40 ju, There could be an 
admixture of cooler grains, To explain the observed 
40 M flux with grains with emissivity similar to Moon 
rock 14321 at 40 m, would require a mass in colder 
grains more than 6 times larger than the dust which is 
responsible for the 10-20 ju flux. 

An alternative explanation is that the observed 
spectrum is contaminated at the longest wavelengths 
by diffuse background radiation from the nearby Orion 
molecular cloud. Maps of this region (Low, Rieke, and 
Armstrong 1973; Werner el al, 1976) indicate that the 
contamination in our beam can be several percent of the 
peak KL flux in the 30-50 n range. To indicate how this 
would affect the observed spectrum, we have subtracted 
3 percent of the KL spectrum (Fig, 1) from the Trape- 
zium data; this is also shown in Figure 2. .After allowing 
for a reasonable amount of molecular cloud contamina- 
tion, the fit between the simple silicate dust model and 
the data is quite good longM'ard of 25 /i. 

The simple model disagrees with the observed spec- 
trum in one further detail: the local minimum in 
emissivity of the lunar silicate around 23 m is not present 
in the Trapezium spectrum, A similar feature is seen 
in all the lunar silicates studied by Knacke and Thom- 
son (1973), Terrestrial silicates also show fine structure 
in the 20 n region. Thus either the fine structure is 
masked by a mbcture of different silicates or the cosmic 
material is different from that in lunar and terrestrial 
silicates. Day (1974) has found a similar lack of struc- 
ti-ire in the 20 n region in an artificially produced 


amorphous silicate and in carbonaceous chondrite 
silicate material, but they peak at longer wavelengths 
than the Trapezium, Recently Penman (1976) has 
measured the optical constants of several meteoritic 
silicate materials; a combination of the Vigarno and 
Murchison carbonaceous chondrite materials would 
peak at 18.5 g and show no further structure in the 
20 ju region, It would seem that these amorphous, dirty 
silicate materials may provide a closer representation 
of the observed Trapezium spectrum, 

V, Tire KLEIN', MA.S’N'»t.OW N'EBUU 

The 16-40 a spectrum of the Kleinmann-Lo\s* Xebula 
(Fig, 1) shows a broad continuum with no sharp 
features, In the region of overlap, the spectral shape is 
in quite good agreement with the data presented in the 
accompanying paper (FS). The approximate equality 
of 20 (1 color temperature and brightness temperature 
with the radio determinations of gas temperatures in 
the molecular cloud (Liszt el at, 1974) leads us to 
conclude that the actual physical temperature of dust 
grains is approximately 80 K in the region of the 
molecular cloud which is seen at these wavelengths 
(FS), The absence of the rotational lines of molecular 
hydrogen near 17 and 28 n is not unexpected, because 
of the coarse resolution (X/AX^ 30) of the present 
experiment. For an optically thick line of width 10 km 
s">, gas at temperatures from 0 to 100 K would result 
in at most a ~0,l percent feature in the present 
spectrum, 

The small residual absorption feature at -^20 m is 
ascribed to the effects of radiation transfer in the cloud 
coupled with the wavelength-dependent opacity of the 
dust grains (FS; Knacke ancf Thomson 1973; Kwan and 
Scoville 1975); because of the reduced opacity at 16 
and 24 p we see to a physically deeper layer in the 
cloud than at 20 n and find that the dust temperature 
is higher there. Most dust models and the spectrum, of 
the Trapezium (Fig. 1) indicate that the dust emissivity 
decreases with wavelength longward of 20 ju< Therefore, 
we should see deeper and deeper into the molecular 
cloud at wavelengths longer than 20 n. The obser\’ed 
fact that the brightness temperature peaks at X — 28 u 
and then decreases to longer w.avelengths indicates 
that at X— ' 28 JU we are seeing to the central, hottest 
portion of the nebula and thereafter we see beyond the 
center where the dust is cooler. Thus we conclude that 
t 5 s» « 1 to the center of the Orion molecular clout*. 

The 28 ji optical depth may be used to estimate the 
total mass and central density in, the Orion molecular 
cloud fOMC-l). Following the discussion of Gillett 
el al. (1975), If the density in the cloud goes as p <x 
1/R* (Larson 1969), the central gas density pn (gas) and 
total cloud mass Mai will be related to the observed 
optical depth t\ by 

and 

- Uf—mx , ( 2 ) 
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where i?» ami «re the Inner and outer radii, 

respectively, of the region obset^-ed, I •/!$ the mass ratio 
of gas to dust, and Hk is the mass absorption of the 
coethcient of the dust. For normal galactic abundances, 
t / ^ 300 for silicate material, For the lunar silicates 
studied by fCnacke and Thompson (1973), na,, m f 
Kio« a* 4 X I0‘cm^ g"‘. In the present case, we identify 
the inner radius with the 12!5 beam radius which, 
at the Orion distance of 500 pc, corresponds to 0.9 X 
10‘^ cm, Then po, e.xpressed as a number density of 
molecular hydrogen atoms, is Ha(Hj) ^ 2 X 10’ cm"’ 
for the central 25" region of the molecular cloud, If the 
outer radius R\ is identified with the -^4' X 9' emission 
ridge seen in the CO and CS molecular maps of Liszt 
it at. (1974), the total mass svithin this region is given 
by Mm > 600 Mq, From considerations of the lor- 
maldehyd^ emission in Orion observed with a 1' beam, 
Harv’ey et at. (1974) estimated an Hi density of ~10’ 
cm“- and a total mass of —500 Mq within a 2' X 5' 
region. The mass estimates are in quite good agreement: 
the larger density inferred here >,ould be due to the 
smaller beam size employed in the present study and 
the assumption here that p « l/R^ in this cloud. The 
assumption of a steep density gradient has recently 
been supported by a study of the 1 mm continuum 
emission from Orion and three other molecular clouds 
by Westbrook el at, (1976). 

Since kio/ksJ), — 7 for lunar silicates, we would also 
infer that no,, — 7 to the center of the molecular cloud, 
Gillett el al, (1975) found a ratio of visual extinction to 
10 M optical depth .l./rio,, > IS from a study of 
compact H u regions embedded in molecular clouds. 


This would imply a visual extinction of >100 mag to 
the central region of the molecular cloud. This amount 
of extinction is in good agreement with that deduced 
from microwave observations of formaldehyde in the 
same region of the Orion Kebula (Thaddeus el ah 1971 1. 

VI. COKCLUSIOKS 

The Trapezium spectrum shows a broad emission 
feature peaking at about 18,5 p. This supports the 
identification of silicates as a major component of the 
dust in this source. At longer wavelengths the spectrum 
may be contaminated by radiation from the nearby 
Orion molecular cloud. 

The Kleinmann*Low nebula shows a smooth spec* 
trum witi: only a small residual absorption feature near 
20 p and a gradual decrease in brightness temperature 
longward oi 28 p. The observed spectrum can be 
understood as a result of radiation transfer in the 
molecular cloud; this indicates the hottest, central 
region of the Orion molecular cloud is being observed 
at approximately 28 p, 

We wish to thank W. J, Moore of the Xaval Research 
Laboratory for providing the GesOa detector which 
gave alwavs welcome increased sensitivity from 20 to 
40 p. 
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rurthtrpsr*, If tht KmiM **C/‘*C r*tlO U In f»8t 
no grenUp thin tho tfctn tht 

prf»«n% rsiul^i p«i *0 yppor *tnit «f *30S on sht 
rr»etian»l *l)t>nfS«n«tt or urgon in tht atMiphtr*, 
Strone tnittion fttturts vtrt found in tht stnttri 
tt tht ?50 e»'* ’*C**0i »h»6pp%ien Jinti, Thttt htvt 
iiUisitn frofilt* *od ieply » ‘glottis ttrjtpttur* of 
mifo K »t Mioiphtrie tltitudtt. Tht nhtolutt 
ietaRtltiit of tht tniitioa ftttuptt wt, hovtvtr, 
frsm tvo to tin tistt srttttr thtn txptettd fron {»• 
in mTE ftt tiiii ttaporttwi, Tii« tsiiiien tppitrt 
to pffult iVo« CO* tMitttion hy loltr ntwf-finfytptd 
p»di*tien. Thii vork vm lupporttd in port hy }(AfA- 
Crunt l!0t-05-003.2t2 wid »OR-05<*<03->‘52. 


32.04.02 Thf rif_|nfrii'td gptct^ o f CuplttP . 

D. 3cc'Rv:tc,h3.f 7 o.V. sJ’&idiSft, sj.mJ30s, 

C.0. 30 AR 3 U, t.c. mom Md y.c. witteiorw, mA-j&<n» 

ytttftroh Ctnttr. Moffttt Titld. CA ?403g. • Tht sptf 
trua of dupittr hot httn Btuurtd fron Co>l20 aicroni 
froa m nititudo of 13 ka hJOOO ft) utlas tht 91»fia 
ttluscopt on hotrd RASA'i Kuiptr Alrhornt Ohtar/ttory. 

Tht «pcetroatt*r it t Hlohtlton InttJrftreatttr vlth t 
uitfttl htndpMS froa Co-350 ca*^. Tht rtsolutlon of 
tht tptetr* is thout 5 ca**. f?oa 50 to 100 siorons, 
tht tuda opacity toupct it dut to the rotttioniil epto- 
trua of }IH, plut the Ha trtntlttiontl csntinuua. 3yn* 
thttie sptetrt btvt httn ctleulatcd uting both tht tffij 
& Ha opacity and tht t(tuatioa for radittivt transport. 
Stvtral rtpontd prttturt itaptrtturt profiltt for tht 
Jovian ataotphtrt and tht saturattd vapor pretsurt at a 
function of ttsptraturt for KK3 art uitd. In. tht 
unsatvurattd rtgions a aiadag ratio of 1.5x10’^ is uttd 
for tht partial prttturt of JIH3. Tor tht'prtssurtt uttd, 
tht iOla lint shape t art given oy the loventsian profile 
Vitb hroadtning duo to K*, Kt and hH* included. She 
synthstis Ipoet*'* msm SSMpsrgd to t ha obigrvatloaa and 
diacucaed. The aaount of MH3 in the linfof-iight it 
deduced. ^ 


32,05.02 Ctnvtt of Orovih for the 0«0 3(J) H» guad- 
rupolt tintt for Jupiter and Orasut. 0. OOCRYIlca and 
C. CKACKStlAh, JR., WA3A-Aati Rtteareh Center, Moffttt 
yield. CA_ 9U035 . - Since high altitude platfoias, tuch 
at the C«i4i, balloont, and the Shuttle art htooBinf 
availahlt for carrying tcletcopet above the tropopeuie, 
the 0*0 S(J) (luadrupola linet of once obscured by 
telluric ^*0 can be observed in the absotphertt of the 
outer planets with hi|p* resolutioo inttruaentt. The 
dcvtlopatp''. iof narrow freciUency lasers operating in the 
infrared hw s Bade posalblt the axtantion of radioaatro- 
nosdcal techniques to this region of the sptetruo. Wt 
have calculated the cum of growth and the line shape 
for the 29, 17, and 12^ Unea for Jupiter and Uranua 
using the Caletry profile. For each line, ve calcu- 
lated the curve of growth using the cool, ncoinal and 
vara models of T.h. Olvlne (HASA SP-8O69) and y.O. 
Palluconl (NASA SP-9103) for Jupiter and Uranua, 
respectively. These three Bodcls repreeent different 


Hj/He ratiei which can be supported by obiirvatlons. 
The ataospherei are di-rtdtd into about 33 l#jcf» 
the abiorptioa it salcuiated for enission up through 
these layers. The line strsneths were taken iFe« the 
Uterature, the diffusiwj constant and erttsure 
breadiivinc coefficient are calculated taking into 
account the intemtiens between hydrogen and heiiua. 
Also a seai»«rFirical lint shift wae included which 
tends to lessen the saturation. We present the curve 
of CTCwth for the strsneeat lint C-o S(l) at I7t» for 
both Jupiter and Uranus noaintl sndel. 
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32.06.02 VtBUSt Tbt 16-10 Mtev on 
REED, CopeU.V._a.aru»w.'. Aircraft. 
tou, J.R, HOUCK a u, - Wt report 

on observations of Venus tstde in the epeetral Intcrvsl 
17-39 Bierone with a resolution of approxlaately 1 Bi- 
cron. The obeervatione were obtained with a two chan- 
nel grtting spectroteter at an altitude of lb ka using 
the RASA Learjet 31 Ca open port telescope. The Venue 
spectrua exhibits three broad and shallow absorption 
features near 17 aierons, 22 Bierone, and 39 alcrons. 
The spectrua le consistent with ealseion froa a hate 
of aqueous sulfuric acid droplets, 


32,07.02 17-25U spectra of Jupiter «nd Ssfcurn. 

A, TCKUIAOA, R. F. XHACXE, S T. OMS, 8 U«-Stony irook- 
Spectra of Jupitac and datum in tha l7-2Su spectral 
region have bean obtained at refolutione of 0.03u and 
O.Sti, tespactively, with a qrou.td-based Hicheison 
speetroaacar. Xheomtlcai spectra ara coopated to 
tha diaarvatlona, Tha shape of tha continuua appears 
to be consistent with recent aCAOtpheric aodali and 
tha preeaure-indueed opacity. of H,. There is no 
conpeilinf evidence for sulfur or other sbsorfaer* 
in Jupiter In this data. 


32. 08.0 2 Far Infrared t.aaellar Gritlnt Observations 
free the Gerard P. Xuiper Airborne ObservecorV . J.I.. 
Fima, M.P. SAVEDOFF, S J.O. DUTHIE, U. Kechester. - 
A Isaellsr irselng incerferoaetat haa been succeesfuUy 
used on the Cerard P. Koiper Airborne Observatory to 
obtain spectra baeveen lOOU and 600U of Jupiter and 
Che Kleinoiann-Low Habula. Va present our first spectra 
of chese.objeeta at a rtsolution of 3.2 ca"! and 
13,1 ca"^ respectively using Mars as a caUbratlon 
standard. The Jupiter spectrua exhibits brightness 
teDperaturas coniiatant with existing broadband oasa- 
urenenta, and we eoapar($ it to thebrecical nodala. 

We failed to detect KH3 dn theea uavclengche. The 
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The Venus emission spectrum was measured from the NASA Tear Jet On five nights in June 
1073. A cooled grating spectrometer with a resolution of X/AA £*23 over the spectral Interval 
17 to 38 unx was used. The main features In the observed spectrum are consistent with the 
theoretical emission spectrum of a haze of aqueous sulfuric acid droplets suspended in a COi 
atmosphere. '* 


INTRODUCTION 

TJbs thermal emiggion gpsctrum of a 
planet containg information on tho chemical 
composition and thermal structure of the 
planetary atmosphere. In the case of Venus, 
the peak surface temperature of ~750®K is 
observed only in the radio region ; the 200- 
260 "K temperatures observed in the ther- 
mal infrared originate from a much higher 
le\'el in the Venus atmosphere. It has been 
proposed (Sill, 1972; Young, 1973, 1974; 
Pollack et al, 1974) that a haze cloud of 
sulfiuic acid droplets exists in the primarily 
COs atmosphere of Venus at the levels 
observed in the thermal u ' ared. The 
present 17-38-um spectrum of Venus will 
allow a further test of this hypothesis. 

OBSERVATIONS 

The spectrum of Venus has been obtained 
from an altitude of 14 km with the 30-cm 

‘ Present address : Gmmnmn Aircraft, Bethpage, 
N. Y. 11714. 


telescope of the NASA Lear Jet. The in- 
strument employed was a cooled grating 
spectrometer utilizing both tho first- and 
second-order spectra from the grating to 
achieve a resolution of 25 over tho wave- 
length range 17-38 /im (Houck et al., 1975; 
Schaack, 1975). The beam size of 3.0 arc- 
min included the whole disk of Venus, The 
instrumental, response was calibrated using 
the emission of Mars in conjunction with 
tho model of Wright (1970), which at this 
time gave brightness temperatures of 250 
and 248“Iv at wavelengths of 20 and 34 pm, 
respectively (Wright, personal communica- 
tion, 1970). The observed spectrum, shown 
in Fig. 1, is the average of five nights of 
observation in June 1975. Data points have 
been taken at a density of three per resolu- 
tion element so that the statistical errors 
may be estimated from the scatter of 
adjacent points. The spectrum is charac- 
terized by an overall continuum level, with 
noticeable absorption shortward of 20 ^m 
and longward of ZOpm, compared to a 
245®K blackbody. This spectral shape has 
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Fig, 1. The comparison of the Venus emission spectrum (points) with the syntlietic HsS04 
spectrum (solid line) ns described in the text. The tr/tangles and circles are from the first and 
second-order spectrometer channels, respectively, 


been confirmed by later observations by 
Forrest and McCarthy, who used a slightly 
modified instrument (Forrest and Mc- 
Carthy, 1970), Mariner 10 broad-band 
measurements give a disk brightness tem- 
perature of 245®K in the 3.5-o5 ;um-band 
(Chase et aU, 1974). The Venera 9 and 10 
spacecraft measured a somewhat lower 
(S-40 Atm) brightness temperature of 232 ®K 
(day side) and 241 ®K (night side) 
(Keldyah, 1977). 

The continuum level observed here im- 
plies a continuous source of opacity in the 
Venus atmosphere over the entire 17-38-jura 
range with increased opacity shortward of 
20 am and longward of 30 Atm. Possible 
sources of opacity at these wavelengths are 
gaseous CO 2 , water vapor, and the cloud 
material. The wings of the COj band 
may provide some of the short wavelength 
opacity but the opacity decreases rapidly 
longward of 17 Atm, Water vapor could 
supply opacity at longer wavelengths but 
the absence of structure, particularly at the 
position of the strong HjO band at ~35.5 
Atm (Traub and Stier, 1976) indicates this 
is not the primary source of continuum 
opacity. In the following section it is showm 
that a haze of sulfuric acid droplets can 
provide the necessary opacity and explain 


the depressions seen longward of 30 ^m and 
shortward of 20 jum. 

SULFUEIC ACID MODEL 

A number of recent observations indicate 
that the Venus clouds are composed of 
aqueous suifimio acid. Among these are; 
(1) the wavelength dependence of polariza- 
tion of visible light (Hansen and Hovenier, 

1974) , (2) the shape of Venus' reflectivity 
in the 3-4 Atm region (Pollack et al., 1974), 
and (3) the S- to 13-Atm spectrum (Sinton 
and Strong, iyflO; Gillett, Low and Stein, 
1908; Hand et <iL, 1968; Samuelson et al,, 

1975) . The far-infrared observations of 
Venus pre.sented here offer another ob- 
servational test of this hypothesis because 
sulfuric acid has absorption ma.xima at 17, 
22, and longward of 30 nm. 

The radiative transfer method of Young 
(1974) was applied to the Venus atmosphere 
in order to generate a synthetic .spectrum. 
This method makes the simplifying assump- 
tions of a constant temperature lapse rate 
and a cloud scale height invariant with 
altitude. The data of Jones (197.5) and of 
Palmer and Williams (197.5) for 75% 
H-SOa were used to obtain the Mie scat- 
tering and absorption cross sections of the 
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luuo (In^plotH. Thirt data «houlcl bo rcpro- 
h'ontative of tho Vomiij clouds because the 
aorosol is beliovccl to bo between 75 anti 
S5% HjSOi (Young;, 1975) and the far 
infrared optical constants of aqueous sul- 
furic acid vary little ttithin this range 
(F, W. Taylor, personal communication, 
1976). It was assumed that tho droplets are 
small, with radius < 1 ;tm, so that scattering 
is negligible at these wavelengths. Such an 
assumption is consistent with the value of 
1 /im for the mean droplet radius as inferred 
by Hansen and Hovenier (1974). Because 
of strong CO 2 absorption near the 15 ^m 
fundamental only observations at 20 nm 
and longer were used to obtain the best-fit 
model parameters. In addition, the model 
was terminated below the T = 2.30°K level 
as this is the freezing point of 75% HsSO^ 
(Young, 1973). The two free parameters in 
the resulting model spectrum are an overall 
normalization constant and Ti, tho tem- 
perature lapse rate per unit cloud scale 
height. It was found that values of Ti 
between approximately 35 and 45®K pro- 
vided the best fit to the observed data. 

The synthetic brightness temperature 
spectrum for Ti » 35“K is compared 
against the data in Fig. 1. It is seen that 
the principal features of the Venus emission 
spectrum can be accounted for by this 
simple model. This supports the sulfuric 
acid hypothesis for the Venus clouds. The 
precise shape of the spectrum, however, is 
sensitive to parameters of the Venus atmo- 
sphere such as the size distribution of the 
aerosol and the vertical variation of the 
aerosol mbdng ratio. It is also sensitive to 
uncertainties in the infrared emission spec- 
trum of Mars. One of us (J.B.P.) is pursuing 
a more detailed analysis of the data. 

Using oiur absolute calibration and the 
spacecraft data for the Venus atmosphere, 
we deduce a pressure level for our spectrum 
of roughly 200 mbar. This figure is con- 
firmed by an alternate analysis using the 
optical properties of sulfuric acid and the 
25- to 75-mbar pressure level for visible 


optical depth unity deterniiiU’d by Hansen 
and Hovenier (1074). (Tho o.'ctinction 
efficiency Q» is between 10 and 15 times 
larger in tho visible than in the 20-40-am 
region.) At this level in the atmosphere, tho 
data of McClatchoy el al. (1975) lead us to 
believe that the 17-20-/<m opacity is due 
to both the 17-/um absorption maximum of 
sulfuric acid and COs. 

The value derived for Ti 40®K de- 
serves further discussion because it gives 
now information about the Venus atmo- 
sphere at deeper levels than are accessible 
at shorter wavelengths. This value is 
approximately t\vice tho value expected on 
the basis of direct extrapolation downwards 
from overlying levels of the atmosphere 
(Lacis, 1975). From this it can be inferred 
that below a level corresponding to P < 200 
mbar, the cloud aerosol may no longer be 
characterized by a simple scale height law. 
The present analysis gives little indication 
of what the actual distribution is, but it is 
interesting to consider that a failure or the 
simple scale height law may be related to 
the detailed cloud structure observed below 
~65 km altitude by the Venera 9 and 10 
spacecraft (Keldysh, 1977). 
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A liquid'helium-cooled grating spectrometer has 
been constructed and used to make low-resolution far- 
infrared spectrometric observations of astronomical 
sources from the U.S. National Aeronautics and Space 
Administration (NASA) Lear Jet. A grating instrument 
was chosen for this application because it is simply 
constructed, is insensitive to vibration, weighs little, 
and is small. These are important considerations for 
observations with the 30-cm Lear Jet telescope. More- 
over, since our detectors behave as though background 
limited, the sensitivity Of a cooled grating instrument is 
competitive with a multiplexing spectrometer. 

The instrument is shown schematically in Figure 1. 
The focus end of the spectrometer is bolted to the cold 
work surhrce of a liquid helium Dewar; the rest of the 
instrument Is covered by a vapor-cooled radiation 
shield. The warmest parts of the frame are kept at a 
temperature < 20 K under operating conditions. This 
keeps the background flux on the detectors, &om in- 
side the instrument, small compared to radiation enter- 
ing through the polyethylene window hrom the atmo- 
sphere and telescope. To eliminate scattered light, the 
inside of the instrument is coated with an absorbing 
paint, and the direct reflection bom the secondary 
mirror to the detectors is blocked by a black mask on 
tlie secondary. 

The spectrometer has a focal ratio of /6.5, and a 40- 
cm focal length. The entrance slit is 2.2 X 2.9 mm^, 
which defines a beam of 4'5 X 6' on the sky when 
used on the Learscope. The focal ratio and length cam 
be conveniently changed by substituting different sec- 
ondary mirrors. 

The instrument frame is constructed entirely from 
aluminum, and occupies a cylindrical volume approx- 
imately 13 cm in diameter and 26 cm in length within 
the Dewar. The primary and secondary pyrex mirrors 
are spherical and coated with aluminium. They are 
spring or clip mounted to avoid breakage during cool- 
down, but ^e instrument retains its alignment well 

•Present Address: David Dunlap (%servatciry, Universiiy of To- 
ronto. Toronto. Ontario, Canada, UC 4Y6. 


even after many cryogenic cycles. The Dewar capacity 
is 2.3 liters, and its hold time is approximately eleven 
hours. 

After a liquid nitrogen precooling phase, three he- 
lium transfers (about 12 liters, including transfer losses) 
over a period of two to three hours are required to 
bring the system to liquid helium temperature. 

The spectrometer is scanned in wavelength by rotat- 
ing the grating with a stepping motor located on the 
outside bottom of the Dewar. The motor connects via 
a fiberglass shaft that runs through a vacuum coupling 
to a pair of bevel gears, a worm, and a sector gear 
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which combine to rotate the grating mount, Tlie spec- 
trai resolution and wavelength coverage are deter- 
mined by the slit sizes and grating blaze angle. We 
have employed gratings with groove spacings ranging 
from 0,10 to 1.73 mm, blazed at ~ 00 jii. Both custom- 
made aluminum gratings and standard Bausch and 
Lomb repUca gratings have been used. The instrument 
has been used over the wavelength range from 45 g to 
113 /i, at resolving powers ranging from 10 to 150. 
The theoretical maximum resolving power of the in- 
strument is ~ 700 at 90 ja, but the limited guiding ac- 
curacy of the Lear telescope has constrained us to use 
fairly large slits and a correspondingly lower resolving 
power. 

The instrument uses two GesCa photoconductive de- 
tectors (Moore and Shenker 1963) which have been 
produced following the procedures of Pipher (1971). 
These detectors are mounted directly to the cold work 
surface of the Dewar, and operate at — 4.2 K. One de- 
tector works in the wavelength range from 45 ja to 80 
H in first and second order of the grating, with a CaFn 
filter; a diamond scatter filter (Armstrong and Low 
1973) cuts off short wavelength light and KRS-5 acts as 
a long-wavelength cut-off filter. The second detector 
works longward of 60 ja in first order, with BaF, plus 
the CaF, and scatter filter blocking light at shorter 
wavelengths. 

The spectrometer employs very simple cooled MOS- 
FET coupled transimpedance preamplifiers (Wyatt, 
Baker, and Frodsham 1974), for which the circuit is 
shown in Figure 2. The MOSFETs and the feedback 
resistors are mountsd as close as possible to the detec- 
tors, Microphonic noise is not a problem in spite of the 
relatively high vibration levels encountered in the air- 
craft. Linearity of the response is sufficient to permit 
use of the moon to give a spectral calibration even for 



Fig. 2 — The simple .MOSFET coupled transimpedance preamplifiers 
used in the spectrometer. 


faint Infrared sources. The accuracy of this approach 
was confirmed by observations of Mars, whose bright- 
ness as a function of wavelength and time of observa- 
tion has been calculated by Wright (1976). The noise 
of the preamplifier is about 1,3 microvolts-Hz"'^-®, 
which is far less than detector noise under normal op- 
erating conditions. 

In laboratory tests without an integrating cavity at 
— 85 fi effective wavelength and backgrounds of ~ 
10~® watts, these detectors gave measured photon 
NEP’s of 1 X 10"'^ w-Hz"”*, This is about a factor of 
10 larger than the theoretical background limit for a 
photoconductor under these conditions. The current re- 
sponsivity was approximately one ampere per watt at 
the optimum bias voltage. 

In actual observations of M42 and Venus on the air- 
craft, the observed system NEP at a resolving power of 
150 was about 4 X 10“‘^ watts-Hz“® ® at 88 /i and 9 X 
10“ watts-Hz"*’ * at 52 ja. No corrections were made 
for telescope, chopper, and atmospheric losses. Noise 
contributions from guiding errors which generally dom- 
inated over the detector noise were included. The de- 
rived system efficiency is ~ 10%, again including tele- 
scope, chopper, and atmospheric losses. 

Normally a calibration spectrum is obtained with a 
laboratory source prior to each Lear Jet flight. For this 
calibration the spectrometer views a liquid-nitrogen- 
cooled blackbody source placed ~ -.70 cm from the 
polyethylene vacuum window.^ The solid angle sub- 
tended by the source has a diameter of ~ 8°. It equals 
tl;e angular diameter subtended by the wobbling sec- 
ondary mirror at the entrance slit of the spectrometer 
diuing flight. Figure 3 compares preflight and in-flight 
spectra obtained on typical days. The inflight spectrum 
was obtained at a cruising altitude of 13.3 km. The 
preflight spectrum was taken at sea level, at Moffett 
Field, California. The spectral resolving power is — 
120 in both spectra. The lunar spectrum shows some- 
what shallower absorption features than the spectrum 
of water vapor in 70 cm of normal air traversed at sea 
level. The water vapor features are pressure broadened 
in the lower ^<tmosphere. A more detailed comparison 
of the two (Rurves would require a curve-of-growth 
analysis of thit water vapor features. 

The spectral resolving power attained by the in- 
strument varies with the source viewed. Compact 
sources whose images do not fill the entire entrance 
slit width produce spectra with a resolving power of 
~ 150 while extended sources like the moon provide a 
resolving power of 120. Similarly a small shift in the 
position of a spectral feature is produced if a small 
source is imaged first at one side and then the other of 
the entrance slit. Both of these problems could be 
overcome to a large extent by using a field lens in 
front of the entrance slit. 
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Flc. 3— Spectrum of the moon as seen' from- a 'Lear Jet altitude of 13.5 km and a preflight calibration spectrum of a laboratory source seen 
through a path length of 70>'cm:of air at sea level. •< 


The reader can gain a clear understanding of the ca- 
pability of the instrument .from published spectra. 
Low-resolution spectra (A/AA - v 10-.30) are given by 
Ward, Gull, and Harwit (1977a,b), High-resolution 
spectral line data are given by Ward et al. (1975) and 
Dain et al, (1978) ((Oni| 88.16 >),■ and Melnick et al. 
(1978) and Melnick, Gull, and Harwit (1978) ([Oiii] 
51.8 n and [O i] 63 /i). 

We would like to thank Martin Harwit, George Gull, 
and Carl Frederick for their extensive contributions to 
the development and use of this instrument. The e.xcel- 
lent support, of the Airborne Sciences Office at NASA 
Ames made the observational phase of this project pos- 
sible. This work was supported by, NASA contracts 
NCR 33-010-146 and NGR 33-010-182. 
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ABSTRACT 

Emission of [S in] 18.71 /xm radiation has been detected in the Orion Nebula, Ml 7, NGC 2024, 
and WSl, and an Upper limit has been set for Sgr A. The line flux from M17 has been mapped 
with a resolution of T.1, A simple theory for predicting the [S lu] flux from radio continuum 
observations is presented and compared with the data. It is found that Orion and M17 are in 
satisfactory agreement with the observations, and that NGC 2024 is consistent with a low excita* 
tion source behind a reasonable amount of 18 /xm dust extinction. In WSl the low observed 
flux may be explained either by iow S'^ abundance or, more likely, by dust extinction. In Sgr A 
attenuation by dust is the most plausible explanation. 

Subject headings: infrared: sources — infrared: spectra — nebulae: general 


I. introduction 

Infrared fine-structure lines can be used to investigate 
the properties of the dust and ionized gas of H n 
regions. If the ionized nebula suffers little visual 
extinction, as in the Orion Nebula, the observed 
fluxes from fine-structure lines can be used to determine 
electron densities, temperatures, and ionic abundances, 
which in turn can be interpreted to yield information 
about the exciting star and elemental abundances. On 
the other hand, in the case of heavily obscured 
nebulae, such as WSl, the attenuRtion due to dust is 
not negli^ble, and the comparison of observed fluxes 
with predictions based on simple models can be useful 
in estimating the optical constants of the dust, which 
at present are poorly known. 

In the mid-infrarra spectra of galactic H il regions, 
the 18.71 /tm line of [S in] is among the strongest 
observed. It arises from transitions between the ’’Pj 
and ^Px levels in the split ground state of the S**"^ ion. 
Collisions with electrons populate these levels, and for 
low enough densities (n, ;$ 6000cm"“) emission of 
18.71 /xm radiation is the dominant de-excitation 
mechanism. For higher densities, de-excitation due to 
collisions with other electrons becomes dominant, and 
the emissivity divided by ionic and electron densities 
drops bff as 1/n,. The line’s location near the center 
of the 16-20 /xm silicate absorption band makes it an 
ideal probe of the properties of silicate grains, thought 
to be a major constituent of interstellar dust. Knowl- 
edge of the strength of the 20 /xm baud would be a 
further constraint on the chemical composition of the 
grains and would allow better theoretical modelf/of 
the spectra of cool stars and H n regions where Hust 
emission is important. Furthermore, <unce the S' f ->■ 
S + + ionization potential (23.3 eV) is close to th: ,t of 
He“->-He+ (24.6 eV), comparison with radiji; re- 
combination line observations can give clues f;J the 
ionization structure. , / 

Earlier work has detected [S m] 18.71 /xm in the 


Orion Nebula (Balutcau et ai. 1976), in W3A (Moor- 
wood et al. 1978), and in G333.6-0.2, NGC 7027, 
and BD+30®3639 (Greenberg, Dyal, and Geballe 
1977), This paper re^rts on detections of [S iii] 18.71 
/xm in Orion, M17, WSl, and NGC 2024, and on an 
unsuccessful search for line emission in Sgr A, We have 
also mapped line and continuum emission in M17. 

U. OBSERVATIONS AND INSTRUMENTATION 

The spectra reported here were obtained over a 3 
week period in 1977 September and October, using the 
NASA Lear Jet 30 cm telescope operating at an 
altitude of 45,000 feet (13.7 km). The circular beam 
was 217 in diameter and the chopper throw was 8(1 
in azimuth, except for the M17 map, where the 
amplitude was increased to 16(2. Except for the Orion 
Nebula and M17, the 8(1 throw was large enough to 
ensure that the reference beam was in empty sky, and 
in those two cases the line flux in the reference beam is 
estimated from radio maps (Schraml and Mezger 
1969) to be 10% or less of the flux in the main beam. 

Spectra of the subsclar point of the Moon were 
used to correct for the instrumental response and the 
transmission of the Earth’s atmosphere by assuming 
that the Moon’s intrinsic spectrum was a blackbody 
of temperature 395 K. The absolute flux calibration 
was based on observations of Mais and Wright’s 
(1976) model for its brightness temperature at 20|xm. 
(A recent more detailed model [Christensen 1979j of 
Mars would lower the calibrated fluxes by 1.17.) The 
flux Calibration was checked by observing VY CMa„ 
the 20 /xm flux of which agre^ with Morrison and 
Simon (1973). The uncertainty in this calibration is 
bard to estimate but may be of order ± 10%. Night- 
to-night gain variations were ±10%, or less. The 
wavelength calibration was based on laboratory 
measurements of water lines and was frequently 
checked by examining the position of zero order. 

During this flight series our limiting source of noise 
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was guiding inaccuracies. Since our dewar samples 
several closely spaced wavelengths simultaneously, it 
was possible to make a correction for this guiding 
noise. This was performed by assuming that when the 
grating was in one position the average of the 10 
channels should be constant within limits set by the 
other sources of noise, e.g., sky noise and background 
fluctuations. When an average fell outside these 
limits, oil contributing channels were multiplied by 
the correction factor necessary to make the average 
agree with a predetermined value. In the case of 
sources at the peak of infrared maps, like WSl and 
M17 S, this value was the maximum observed, while 
for sources like M17N and Ori, which are on the 
slopes of much stronger sources, the averages were 
corrected to agree with the mean of all the averages. 
Typically these guiding corrections were made to 
about one*ihird of the integrations and were of order 
IQ^^g. They imi)roved the shape of the spectrum, thus 
aiding in Identifying the line. However, these correct 
tions also introduced an additional uncertainty into 
the absolute flux calibration. 

These observations were carried out with a new 
lO^hannei liquid-helium-cooled grating spectrom- 
eter. The 10 Si:Sb photoconductor detectors are 
mounted in the focal plane of the f/7 single pass 
Ebert-Fastie configuration, and together they cover a 
S nm band with a FWHM spectrm resolution of 0.5 
fjim. The entire 16-30 region can be sampled at a 
density of three points per resolution element with a 
total of only nine grating movements. From observa- 
tions of VY CMa, the in-flight system NEP at 23 /xra 
has been estimated to be 8 x 10“^* W in each 
0.5 nm bandpass, with all the detectors equaling this 
performance within ± 10%. By allowing for spectrom- 
eter transmission, beam splitter reflectivity, and 
chopper efSciency, we estimate the in-flight detector 
NEP to be 1.3 x 10"^* W Hz"^'® for each detector. 

Detector output currents first pass through G118 
FET transimpedance amplifiers on the cold spectrom- 
eter baseplate, thereby avoiding long high impedance 
leads which are susceptible to microphonic pickup. 
The signals are next amplified by operational amplifier 
(op amp) circuits with several available gain settings 
and then filtered by Jew pnss filters with 100 Hz 
roll-off. At this point the 10 channels are also averaged 
by an op amp summing circuit, and this “ sum channel ” 
is synchronously demodulated by a lock-in amplifier 
and displayed on a strip chart. In the absence of 
correlate noise, this channel gains a -^10 improve- 
ment in signal-to-noise ratio and can lx very useful 
in finding source peaks. All data channels are fed into 
an IM9A1 8080 microcomputer where they are 
sampled by a 16-channel 12 bit analog to digital 
converter. The synchronous demodulation process is 
carried out in software, using an INTEL 8080 
assembly language program, which normally samples 
and demodulates 12 channels in 1.56 ms (640 Hz). The 
phase lags between the chopper reference signal and 
the detector outputs are the same within ± 3°, and the 
synchronous demodulation routine will accept f», 
sintfie phase delay of up to 32 ms or 180°, whichever 


is smaller. The computer can also step the grating 
according to an easily modified table in its memory, 
and the results of all integrations are written onto 
floppy disks and printed on a small data logger. The 
observer controls the system by means of sense 
switches which call master routines, and by using a 
keyboard to enter alphanumeric data. At all times the 
program's status is displayed on a television monitor 
and on the front panel LEDs. On board the Kuiper 
Airborne Obs°irvatory, the entire data-acquisition 
process is automated, with the IMSAI acting in the 
role of master computer and using the on-board 
minicomputers as slaves for real-time data analysis 
and for nodding the telescope. The data-acquisition 
propam occupies 8K bytes of memory plus an 
additional 4K bytes for a disk operating system. 

UI, RESULTS 
a) Full Spectra 

[S lu] 18.71 /im emission was detected from three 
beam positions in the Orion Nebula (Trapezium, I^, 
and Ori), M17 S, M17 N, NGC 2024, and W51. 
No line emission was seen in Sgr A. Between 17.2 and 
20.2 pm, the spectra were fitted to a quadratic base- 
line and a Gaussian line profile of FWHM - 0.5 pm 
and central wavelength Ag > 18.71 pm. The results are 
plotted in Figure 1 and summarized in Table 1. Where 
large enough to plot, the errors shown in individual 
data points are standard deviadons of the mean 
derived from the scatter of separate integrations. The 
quoted errors on line fluxes derive from the fitting 
procedure and do not include uncertainties in calibra- 
tion or gain changes. Tlie relatively high x® values in 
the Orion and M17 S fits can be attributed to two 
causes. First, there appear to be small systematic 
point-to-point variations in the spectra which show 
up as jitter and scatter and which are not described by 
the statistical errors. Such scatter might be caused by 
incomplete cancellation of the instrumental response 
or by guiding noise. If the statistical error of the 
spectnJm of the Trapezium were increased to 2% of 
the flux, the x® of the fit would be decreased to a 
more acceptable 1,05. Second, in these high signal-to- 
noise spectra the line profile suggests a non-Gaussian 
shape. No attempt was made to eliminate these 
effects, since it was felt that the desired quantity, the 
line flux, was sufficiently well determined. 

b) M17 Map 

Fourteen positions in M17 were examined for [S ni] 
emission; the line was detected in 12 of them. The 
positions were separated by 118 and are roughly 
Oriented along and perpendicular to the axis between 
M17 S and M17 N. Offsetting was from the star 
BD - 16'’4816 by means of an xy stage. The chopper 
throw during the mapping was 1612, ensuring that the 
reference beam was completely off the H n region as 
shown in the radio maps (Schraml and Mezger 1969). 
The sp^ra were sampled at a density of one point per 
resolution element, and due to the lower density the 
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iitting function was simplified to a linear baseline pius 
a Gaussian line profile with FWHM ■» 0.5 {xm and a 
central wavelength of 18.71 /tra. The results arc 
listed in Table 2 and plotted in Figures 2b and 2c, 

IV, CALCULATION OF EXPECTED LINE STRENGTHS 

In the absence of extinction, the expression for the 
[S III] flux from a cloud of Ionized gas may be written 
(Simpson 1975) 

n, 

where is the cmissivity of the ion, 

n„ n, are the number densities of the S** ions, 
electrons, and protons, £ is the emission measure of 
the region integrated over the beam, and 0. Is the 
solid angle of the beam. 

The emissivity of the ionized gas depends on the 
electron density and temperature. For ;$ 6300 cm“®, 
collisionai de-excitation Is negligible, and the emissivity 
is essentially independent of density and only wealciy 
dependent on temperature (Simpson 1975). Due to 
collisionai de-excitation, the emissivity decreases with 
increasing n, for n, > 6300 era'® and falls to hulf of 
itUow density value at n, » 2.5 x 10* cm"®, Observa- 
tions of Sp A fekers ei ai. 1975; Pauls et aJ, 1976) 
indicate that coUisional de-excitation is negligible In 
that source, According to Schraml and Mezger 
(1969), the electron densities of all our Other regions 
are also below 6300 cm"®, However, in at least one 
case, W51, higher-resolution observations show that 
most of the radio continuum flux originates from 
compact, high density sources, On the other hand, 
high density sources are ebsent or a small fraction of 
the total flux in M17 (Matthews, Horten, and Goss 
1978; Webster, Altenhoff, and Wink 1971) and in’ 
Orion (Webster and Altenhoff 1970), The effects of 
high density regions on individual sources are con- 
sidered in more detail below. For the purposes of our 
simple theory we have assumed T, — 10* K and n„ ;$ 
63C0cm“® throughout, which gives ^ 

5.2 X 10"®® ergs cm® s"^ sr"^ 

The abundance of S+ + was taken to be the average 
of 12 positions in the Orion Nebula which were 
measured by Pcimbert and Torres-Peimbert (1977) 
giving n^jxJnp » 1.4 X 10"*. When corrections for 
unobserved ionization states (e.g., S rv) are made, it 
is found that roughly half the sulfur is in S'*' most 
of the rest being S *"*"*■. The variations within the 
nebula are roughly a factor of 2 above and below the 
mean. For the purposes of this simple theory, we 


Fio. 1, — Observed [S m] 18.71 M>n lines. Sampling Is three 
data points per resolution element, and AA » 0.5 >im. The 
error bars are 1 standard deviation of the mean, and the 
solid lines are fits to the dau. The spectra have been normalized 
for legibility. Object and normalization factor: (a) Trapezium 
X 1.33, (6) KL Nebula X l.(c) F* Ononis x 2.(d)M17S x 1. 
(«) M17 N X 1, (/) W51 X 1,123. (^) NGC 2024 x 0.9, (A) 
Sgr A X 0,33. The curve through the Sgr A data represents 
the fitted value of the line strength, which is much less than a 
3 or detection. 
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Fio. 2. — (d) Featufcs in the H n region M17. Aattrisk, the foreground star BD 16*4816 which was used for otfsei guiding. 
Squart and opn trianglt, positions of radio peaks G15, 0-0.7 and GlS.l -0.7 (Scbraml and Mezger 1969). SoUd lint, haif>power 
contour of 15,4 GHz map, 2' beam (Schraml and Mezger 1969), DoUtd lint, half-power contour of emission (Lada, 

Dickinson, and Penfieid 1974). fllltd triangli, peak of Ha emission (Dickei 1968). Hatehtd elllpsi, center of heavily obscured star 
cluster (Bectz tl al. 1976). Numbers give values of A* to H ii region as derived by Dickei (1968) by comparison of tadio and Ha 
fluxes. (A) Map of the integrated [S m) 18.71 /nn line fltu. Units are 10~^* W cm"% and typical errors are ± 1.3 x 10"^* W cm~*. 
(c) Map of the 18.71 pm continuum. Circles represent the beam positions, and the central numbers give the flux in the 2(7 beam at 
that position. Units of flux are W cm*’ Mm*S and typical errors (1 standard deviation, derived from the fltting procedure) 
are ±0.7 x lO"** W cm*’ (d) Map of the ratio of the observed [S nij 18.71 pm line flux to that predicted by the theory 
described in the text. Typical errors in the ratio are ±0,12, 

havb assumed that this abundance is constant for ail 
the regions observed. Finally, the ratio rttln, was taken 
to be 1.10, since most of the helium should be singiy 
ionized in the region where S+ + is found. 

Except for Sgr A, the emission measure was 
estimated from the maps of Schraml and Mezger 
(1969). By assuming that the brightness temperature 
was constant over the Lear Jet beam, the antenna 
temperature at the cc.<ter of the beam was used to 
calculate the radl r Sux in the Lear Jet beam. The 
resulting fluxes are given in Tables 1 and 2. This 
method overestimates the flux from a point source and 
underestimates for a source about the size of the radio 
beam (2'), The error is believed to be no more than a 
factor of L5. 

In the case of Sgr A, aperture synthesis observations 
(Ekers et d. 197S) have separated the thermal source, 

Sgr A West, from the probable supernova remnant, 

Sgr A East. The integrated flux from Sgr A West is 
26 ± 4 Jy at 5 GHz, and since the FWHM of the 
source is only I' x <0^6, ail of this flux was included 
in our Lear Jet beam. 


Given the radio flux in the beam, the emission 
measure can be calculated if a value for T, is assumed. 
Under these assumptions, one can then write a simple 
relation between the [S lii] 18.71 pm flux and the 
optically thin radio free-free continuum, 

,F(18.71 pm) - 2.04 x 10-»/(r„ v)(S,/l Jy)Wcm-» 

where 

/(r„ y) - (r,/10* K)‘'“ an (4.96 x 10^r.=’'»v-0]"' , 

and is a slowly varying function of the electron 
temperature and the radio frequency. For T, ■■ 10* K 
and V- 10 GHz, /(r„ v) - 0.118. The total un- 
certainty due to temi^rature is small. In the low density 
limit Tt - 5000 K would lower the predicted flux by 
5%f while r, - 2 X 10* K would nuse it by 10%. 

Several of our sources lie behind large column 
densities of dust, and if this dust contains silicates one 
can expect attenuation of the [S ill] radiation by the 
18 pm resonance. Because the strength of the 18 pm 
feature is uncertain, we have chosen not to include this 



No. 3, 1979 


£S nd 18,71 MICRON EMISSION 


717 . 



extinction in our colcuk,tions but instead to consider 
the 18 urn optical depth as a quantity which can be 
measured by companng the observed flux with that 
expected in the absence of extinction. 

Hie results of these calculations are given in Tables 
1 and 2 and graphed in Figures 2d and 3. 


V. DISCUSSION 

There is a systematic trend for the predicted 
[S in] fluxes to be too high by an average factor of 
about 1.3, excluding W51 and observed upper limits. 
This could be due to one or more of several effects. 
First, our flux calibration could be low, but the 
agreement of our spectrum of VY CMa with ^ound- 
based measurements argues against this being the 
major factor. Second, the method of estimating the 
radio flux may have systematic errors. In order to 
estimate the magnitude of such errors, models of 
sources with Gfaussian surface-brightness profiles were 
constructed and the computed fluxes (estimated using 
the method outlined above) were compared to the 
exact values. For the source and beam sizes involved, 
the estimated flux was low by 12% at most, which has 
the sense of increasing the disagreement. Finally, the 
assumed ionic abundance may be too high by a factor 
of 1.3. Given the variations of a factor of 2 in the 
Orion Nebula (Peimbert and Torres-Peimbert 1977), 
this seems quite possible. 

Orion.-— Since the 20 silicate feature is seen in 
emission from the Traperium (Forrest, Houck, and 
Reed 1976), the dust optical depth is small, and these 
observations serve as a useful check on the consistency 
of the theory. In all three of our beam positions the 


agreement between predicted and observed fluxes is 
good. 

0\ia flux from the Trapezium corresponds to a ' 
surface brightness of (4.9 ± 0.2) x 10“* Wcm"®sr“^ 
averaged over our 2(7 beam. TTiis value is a factor of 
1,75 times higher than observed by Baluteau et at. 
(1976) (see also Moorwood et al. 1978) with a 55* 
beam and does not require a lower value of the sulfbr 
abundance than assumed by Simpson (1975), For a 
2(7 circular beam centered on the Trapezium, Simpson’s 
model predicts a [S in] flux of 27.4 x 10”“ W cm”® 
(Simpson, private communication), which Is also in 
good agreement with oiir observations. 

this source the agreement between 
predicted and observed fluxes is good, The [S iii] flux 
generally follows the radio flux fairly closely with the 
possible exception of position 5, Here the beam is at 
the interface between the H li region and the molecular 
cloud, so the excitation may be lower or the extraction 
may be greater. At most of the other points the pre- 
dicted flux is within 25% of that expected, and the 
variations are within the errors. There is no strong 
evidence for systematic changes in the ratio of the 
predicted to the expected (S ni] flux over the nebula. 

NGC 2024.— This source appears to have low 
excitation. The radio recombination lines of He are 
anomalously weak and the size of the He* source Is 
significantly smaller than that of radio continuum # 
source, as well as being offset by about 1' (MacLeod, 
Doherty, and Higgs 1975). Optical observations 
(Balick 1976) show that forbidden lines of high 
ionization spwies are weak or absent, and that the ‘ 
ratio of He* /He integrated over the nebuia is less 
than 0.2. This evidence indicates that the source of 
the ionizing radiation is cool. Balick models the 
H n region with n, »• 10® cm”® and a single star of 
r - 33,500 K. 

Our predicted flux is approximately twice what is 
observed, but the interpretation is conAised by the 
presence of substantial dust column densities in the 
line of sight, Grasdaien (1974) measured "• 8,3 ± 

0,2 mag toward NGC 2024 No. 1 and «- 32 ± 6 
toward NGC 2024 No, 2. NGC 2024 No. 1 is on the 
edge of the dark lane overlying the core of the H n 
region, so its visual extinction may be used as a lower 
limit to the mean extinction of the source. We can 
estimate an upper limit to the 18.71 fxm optical depth 
by assuming that all the sulfhr is in S** and that the 
total sulfur abundance is S/H - 2.6 x 10"* flog S >■ 
7.41 by number] as given by the average of the posi- 
tions in Orion (Peimbert and Torres-Peimbert 1977). 
From this we find t (18.71 ^ra) ;S 1.3, which leads to 
AvIt{IS. 71 fim) ^ 6. The ratio AJr(9,7 fitn) toward 
VI Cyg No. 12 is 14 (Gillctt et al. 19756), and for most 
silicates the 2C ,am resonance is weaker than the 10 jum 
one, so this lower limit on .4,/t(18.71 /xm) is entirely 
reasonable. 

Using Balick and Sneden (1976) to derive S**/S 
based on Balick’s (1976) model of NGC 2024, one 
finds that S**/S»0.7 and T<18.71/xm)» 1. The 
observed [S m] flux is therefore consistent with the 
Balick model. 
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For this giant H ii region, our simpSe theory 
predicts 5 times more [S m] Hux than observed. At 
least three effects can be proposed as the cause of this 
discrepancy: coilisional de>excitation in compact 
regions, low abundance, and attenuation by 
dust, 

Much of the radio flux in our beam probably comes 
from the two compact components seen in the high- 
resolution radio observations of Martin (1972), Felli, 
Tofani, and D'Addario (1974), and Scott (1978). 
Thr.?.? compact H u regions will be referred to by the 
nat'ies given to them by Wynn-Williams, Becklin, and 
Neugebauer (1974), W51 IRS 1 and W51 IRS 2. In 
order to estimate the [S m] 18.71 fim line fluxes 
expected from these sources, the observed radio 
fluxes have been fitted to a thermal bremsstrahlung 
model to derive the “optically thin" flux at 5 GHz. 
This leads to an estimate of approximately 45 ± 10 Jy 
from W51 IRS 1 and 14 ± 2 Jy for W51 IRS 2. 
Scott (1978) estimates the electron densities to be 
n, « 10* cm““ for W51 IRS 1 and n, « 2.5 x 10* 
cm**’’ for W51 IRS 2. Then, using the formulation 
of Petrosian (1970) and Simpson (1975) for the 
expected stren^h of the [S iii] 18.71 Mm line with a 

5 Ui abundance of 1.4 x lO"* with respect to protons, 
the [S in] 18.71 Min line flux should be approximately 
8.9 X 10-i*Wcm-» from W51 IRS 1 and 1.7 x 
10-** W cm-» from W51 IRS 2 or a total of ~ 10.6 x 
10"^*Wsm-® for the compact components alone. 
The observed line flux of 3.3 ± 0.6 x 10"^* W cm“® 
is approximately a factor of 3 lower than this. There- 
fore density effects alone are not enough to explain 
the smaller observed [S ni] flux. 

Due to the enormous visual obscuration of W51, it 
is not possible to directly measure the S*"** abundance 
as in Orion. Instead one must rely on indirect methods. 
Based on continuum measuy ments, Beiging (19)‘5) 
finds that one or more 04 stars are the exciting stars 
in these sources. 04 V stars have temperatures of 
about 50,000 K (Panagia 1973), which is hot enough 
that about 90% of the sulfur is in (Baltck and 
Sneden 1976). One can next try to observe (S iv] 10.5 
Mm to check on this. Using Simpson's (1975) emissivity 
for T, » 10* K and n, » 10'' cra““, a total sulfur 
abundance of 12 + log N, 7.4 (Peimbert and 
Torres-Peimbert 1977), and Martin’s emission measure 
of 72,5 X 10* pc cm-« and size (14' x 31') for IRS 1, 
we predict F (£S iv] 10.5 Mm) =» 3.4 x 10"^®/e"’, 
where is the attenuation due to dust and /is the 
fraction of sulfur in S If we take /- 1 and 
T » 3.9 (Gillett et al. 1975a), we predict F ([S r/] 
10.5 Mm) - 6.8 x lO-^^ W cm"®. GiUett eta/. (1975a) 
did not detect [S rv], and an upper limit of about 

6 X 10“‘®Wcm"“ can be set to the line flux. It 
appears that S is less abundant th,m exacted. 
This is puzzling because [O m] 88 Mm is observed at 
close to the expected strength, assuming an O*"*' 
abundance like Orion’s (Dain et ai. 1978), but O'*"*’ 
and S’*"*"*’ require similar photon energies for their 
creation. Either S v or S in must be the dominant form. 
Orion is ionized by 9* C Ori, which is an 06 or 05 
star (Bedv^f 1964; Smithsonian Astrophysical Ob- 


servatory 1966), so according to Baiick and Sneden, 
> 0.8 and S++/5 ;$ 0.05. Yet the (S nij flux 
from Orion requires that much of the S be S ■*"'■, and 
the [S iv] 10.5 Mm line is undetected and weaker than 
predicted by a factor of 4 (Moorwood et al. 1978). By 
analogy it is plausible that S*+ is the dominant ion 
of sulfur in W51 as well. It is uildear why this should 
be so. Possibilities include (1) the stellar atmosphere 
models are incorrect; (2) selective ultraviolet extinc- 
tion by dust cools the radiation field in the H ii 
region; (3) charge exchange reactions are important 
for sulfur but not for oxygen; and (4) W51 is ionized 
by a cluster of stars of spectral type later than 04. 
For a discussion of this point, see Rank et al, (1978). 

The visual extinction toward W51 IRS 1 and IRS 2 
is estimated to be more than 60 and 24 mag, respec- 
tively, and both sources show deep 10 Mm silicate 
absorption features (Gillett et al, 1975a), Our con- 
tinuum spectrum Shows a shallow 20 Mm absorption 
feature, so large column densities of dust are clearly 
present. If we assume that the observed flux is io’Ver 
than expected because of dust extinction, we derive 
t(18,71 Mm) “ 1 or greater, over most of the source. 

Sagittarius A . — In this H H region, which is believed 
to lie at the galactic center, [S iii] emission was not 
detected and a 3 <r upper limit of 2. 3 x 10-'*Wcm“® 
may be set, which is a factor of 2.5 below the expected 
flux. From radio synthesis maps of the thermal 
component, Sgr A (West), Ekers et al. (1975) derived 
n, - 10® cm"®, well below the value of 6000 cm-® at 
which coilisional de-excitation begins to lower the 
emissivity. Plausible arguments that most of the sulfur 
is in S'* *' may be made from other observations. 
Mezger and Smith (1976) found that the integrated 
strength of the He 109a recombination line was 
0.095 ± 0.04 times that of H, indicating that almost 
all He is singly ionized as long as the He abundance 
is normal. Similarly, observations of [Ne ii] 12.8 Mm 
(Aitken et al. 1976; Wollman et al, j 977) are consistent 
with most Ne being in Ne'*’, unCer the assumptions 
that the Ne abundance is normal and that r(12.8 Mm) 
« 1. On the other hand, [O ui] 88 Mm radiation is 
weaker than expected, indicating that O'*"*’ is under- 
abundant relative to Orion (Dain et al 1978). On the 
basis of the coincidence of ionization potentials, one 
can argue that the radiation field is hot enough to 
make Nc'*’, S’*’ ’*’, and He’*’ the dominant ions of their 
elements, but that the ultraviolet flux is not bard 
enough to create abundant O’*"*’ or to turn S’*"*’ into 
S’*"*"*’. Therefore we feel our as.<Umption that most of 
the sulfur is in S’*’ ’*’ is justified. By elimination, one is 
left with attenuation by dust for an explanation. The 
16-20 Mm continuum absorption feature is the de«ipest 
one ever seen in our Lear Jet observations, leading one 
to suspect that extinction is important here. If we 
interpret the discrepancy between predicted flux and 
observed upper limit as due solely to dust extinction, 
then we derive t(18.71 Mm) > 0.9. Further observa- 
tions with a smaller beam (McCarthy, Forrest, and 
Houck 1979; allow setting a more stringent limit 
to the [S in] flux and permit comparison with other 
wavelengths. 
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VI, CONCLUSIONS 

[S III] 18.71 urn radiation has been detected in the 
Orion Nebula, MI7, NGC 2024, and W51, but not in 
Sgr A. A simple method of predicting the [S iii] 
flux from the radio continuum flux has been applied, 
and it has been shown that Orion and M17 agree with 
the expected flux. The observed line flux in NGC 
2024 is consistent with Balick's model for the exciting 
source, plus a moderate amount of extinction by 
silicate dust grains. In the case of W51, the low 
observed flux may be due to either a lower S+ + 
abundance than in Orion or dust extinction or both. 


The most plausible explanation for the low [S in] 
emission from the galactic center Sgr A is attenuation 
by dust. 
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ABSTRACT 

Airborne observations of the 16-39 /xm spectra of 10 oxygen-rich stars with excess emission in 
the infrared have been obtained. The stars show excess emission attributed to circumsteliar dust 
grains in the 16-39 /xm region in the form of a broad hump peaking near 18 and falling 
smoothly to longer waveien^hs. The emission is similar in character to the emission from the 
Trapezium region of the Orion Nebula, indicating the grain materials are quite similar in these 
objects. ‘l"he existence of a feature in the 20 /xm region is consistent with the O-Si-0 bending 
resonance existed for silicate material. The lack of any sharp structure in the spectra indicates 
the silicate is in an amorphous, disordered form. A simple model of small grains of carbonaceous 
chondrite siliqate material in a diffuse circumsteliar enveloi^e is shown to give a good qualitative 
fit to the observed 8-39 /xm circumsteliar spectra. Comparison of the observed spectra with the 
model spectra indicates the grain emissivity fails as 1/A^ from 20 to 40 /xm. 

Subject headings: infrared: spectra — stars: circumsteliar shells — stars: late-type 


I. INTRODUCTION 

At wavelengths longward of 8 /xm, the spectra of 
many cool stars show an excess in flux which has been 
attributed to circumsteliar dust emission (Giilett, Low, 
and Stein 1968; Woolf and Ney 1969). The circum- 
stellar nature of this excess emission has been con- 
firmed by spatially resolving the shells from the star 
(Zapalla et al. 1974; McCarthy, Low, and Howell 
1977; Sutton et al. 1977) and its identification with 
dust supported by the lack of any fine structure in 
high-resolution spectra (GanunoJ, Gaustad, imd 
Treffers 1972; Treffers and Cohen, 1974). Less certain 
are the chepiical and physical composition of the dust 
gicains, their temperatures and spatial distribution 
about the star, and the wavelength dependence of the 
opacity of the grains. These characteristics are im- 
portant to an understanding of the process of grain 
formation and mass loss from the stars. Further, 
infrared spectroscopy has shown that the grains 
around oxygen-rich stars are quite similar to an 
important component of the interstellar dust grains. 
Thus information on the composition and opacity of 
the circumsteliar dust ^ains has important applica- 
tions to the understanding of other galactic infrared 
sources. 

The present study of the 16-39 /xm spectra of 
oxygen-rich stars provides new information about the 
dust grains in the cimcumstellar envelopes. If the 
10 /xm emission feature in circumsteliar spectra is due 
to a Si-0 stretching resonance in silicate grains, as 
suggested by Woolf and Ney (1969), a second O-Si-0 
bending resonance is expected in the 20 /xm region. 
The shape of the 20 /xm resonance depends on the 
chemical/physical composition of the silicate material. 
For instance, well-ordered silicates, as are found on 
the' Earth and Moon, all show fine structure in the 


20 /xm region. The instrumentation and observations 
are presented in §§ II and III and discussed in § IV. A 
simple shell model of Silicate grains is compart with 
the observed spectra in § V. 

II, INSTRUMENTATION 

The observations were made with helium-cooled 
Ebert-Fastie spectrometers, using telescopes iboard 
NASA aircraft operating between 41,000 and 45,000 
foot (12,505 m and 13,725 m) altitude. The twq- 
channel spectrometer covers the wavelength range 
16-39 /xm (Forrest, Houck, and Reed 1976). A Si: As 
detector was used for the 16-23 /xm band with a 
resolution 0.5 /xm FWHM. Prior to 1977 the 21-39 /xm 
band used a Ge:Ga photoconductor built by W. J. 
Moore with a resolution of 1.2 /xm FWHM. In 1977 
February this was replaced by a more sensitive Ge:Zn 
detector configured to give a resolution of 1,0 /xm 
FWHM. The 10-channel spectrometer covers the 
wavelength range 16-30 /xm and employs 10 Si:Sb 
photoconductors with an Intel 8080 microprocessor- 
based data system which demodulates the sipals of 10 
detectors simultaneously. This spectrometer and the 
data system are described more fully in a separate 
paper (McCarthy, Forrest, and Houck 1979). The 
spectral resolution was 0.2 /xm FWHM for all the 
observations e,\cept for the spectrum of VY CMa, for 
which AA » 0.5 /xm FWHM. The spectrometers were 
used in conjunction with the 91 cm telescope of the 
NASA Gerard P. Kuiper Airborne Observatory 
(KAO) and the 30 cm telescope of the NASA Lear Jet 
(VY CMa only). The beam size of approximately 30* 
on the KAO and 2'.7-3'.2 on the Lear Jet was large 
enough to include most of the flux from these compact 
sources. Standard sky chopping using the oscillating 
secondary was employed. 
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TABLE 1 


Stars Included in the Observing Program 



Name 

(1) 


IRC 

(2) 

AFGL 

(3) 

Spectrum 

(4) 

Variability 

(5) 

Bolometric 
Excess * 
(6) 

® Ori# m • • f r • ( ii » t 

* ( » * ■ 



+ 10100 

836 

M2 lab 

SRc 

s 

X Her 



•*•*«**># 

+ 50248 

1 » # 

M6e 

SRb 

s 

/i Ocp« »#*»»#••• 




+ 60325 

2802 

M2ela 

SRc 

s 

UOO» #)#•>#> 



***#***•« 

+ 30257 

1706 

M7e-M8e 

SRb 

s 

RCas.... 

* t 1 • 



+ 50484 

3188 

M6e-M8e 

M 

s 

NML Tall 

* * * * 

• X • * » • • 


+ 10050 

529 

M6e-M10e 

M 

m 

PZ CTas • » • . t ( « » « 

* • • < 

• X f . X . 

»••»««»#« 

+ 60417 

3138 

M3 la® 

SRa 

m 

VYCMa., 




-30087 

1111 

M5 Ib pec 

Lc 

1 

IRC +10420.,., 

v # • • 



+ 10420 

2390 

F8-G0I" 


1 

OH 26.5+0.6,., 





2205 

... 


1 


* Tlie bolometric excess is the fractional amount of the total luminosity which appears in excess of that expected from the star 
alone as described in the text; s i> smalt ■ 0 to 0.05 ; m « moderate « 0.05 to 0.5; 1 » large 0.5 to 1 .0. 

Note. — The spectral type and variability type are taken from Kukarkin tt al, (1969) except as noted: * Humphreys 1970; 
" Humphreys tt al. 1973, 


in. OBSERVATIONS 

The progratT stars are listed in Table 1 along with 
their 2 micron sky survey number (IRC, Neugebauer 
and Leighton 1969), their AFGL sky survey number 
(AFQL, Price and Walker 1976), the spiral and 
variable type, and an estimate of the fractional amount 
of the total emission which is in excess of that expected 
from the star alone (“bolometric excess”). The stars 
were selected on the basis of their brightness at 20 /xm, 
the presence of the 10 /xm “silicate” feature, and 
spectral evidence in the visible of an oxygen-rich 
atmosphere (0/C > 1). With the exception of IRC 
+10420 all the stars are luniinous, red-giant stars 
probably in an advanced stage of stellar evolution. 
IRC +10420 is an unusual and luminous GO I star with 
large excess emission from dust (Humphreys et al. 
1973). OH 26.5+0.6 shows evidence of having the 
thickest dust shell of the stars in this group, with 
strong^ absorption at 10 /xm and a very cool energy 
distribiition. It is thought to be a Mira variable red 
giant which is losing mass rapidly (Forrest et al. 1978). 
The stars are arranged here in order of increasing 
optical depth of the circumsteilar material, as judged 
from their fractional bolometric excesses (column [6]). 
The method used in deriving this “bolometric excess” 
is described more fully in § IV. . 

The observations reported here were obtained 
between 1976 January and 1978 May. The line-of-sight 
water-vapor column density was monitored aboard the 
KAO and ranged from 5 to 20 /xm precipitable water 
vapor. The absence of water-vapor features in the 
derived spectra shows that the absorption features 
have been adequately removed by the normalization 
procedure. Shortv'ard of 17 /xm, absorption due to 
terrestrial COj affects the spectrum, and it was neces- 
sary to make small differential corrections to account 
for this. For spectral calibration, the Moon was 
observed with the same instruments. The Moon 
was assumed to radiate as a blackbody at a tempera- 
ture appropriate for the phase of and position of the 


beam on the Moon as judged by the 10 /xm measure- 
ments of Geoffrion, Komer, and Sinton (1960). For 
the absolute flux calibration on the KAO, ground- 
based measurement of the 20 /xm flux from non- 
variable sources were used as follows: 1976 May, NGC 
7027 (Fao«m * 6.2 x 10“^® W cm“* /xm“S Becklin, 
Neugebauer, and W|ynn-Williams 1973); 1976 Novem- 
ber, a Ori (Fao,,m =* 1-43 X 10"^® W cji““ /xm“^, 
Morrison and Simon 1973) ; 1977 June, /x Cep (Fao «n “ 
5.8 X 10~“ Wcm““/xm"S Morrison and Simon 
1973); 1978 May, a Her (Fao«m - 3.7 x 10-‘« W 
cm"® /xm■^ Morrison and Simon 1973). For the Lear 
Jet and 1978 January KAO observations. Mars was 
used as the flux standard at 20 /xm, using the model of 
Wright (1976). 

The spectra derived under these assumptions are 
shown in Figures la-le. The flux levels have been 
normalized for convenient display. In Table 2 the 
spectra are listed. Column (1) gives the number of the 
spectrum which identifies the spectrum in the appro- 
priate figure. Columns (2) and (3) give the name of the 
star and the date of the observation. Column (4) gives 
the spectrometer and telescope used, and column (5), 
labeled Fao »n» gives the flux level of the horizontal tick 
next to the identification numbers in the figures; it is 
also the approximate flux level observed at 20 /xm. In 
the cases where the infirared excess is small (Table 1), a 
considerable fraction of the flux observed is from the 
stellar photosphere. To estimate the amount expected 
from the stellar photosphere alone, photometry from 
0.84 lo 5 /xm was fitted with blackbody of tempera- 
ture appropriate to the spectr al type and extrapolated 
to 20 /xm. This flux will be , denoted Fao«m(star); the 
ratio F 2 o«m(star)/F 2 o«m «s giimn in column (/5). As an 
example of the probable contribution of the stellar 
photosphere, the blackbodies for RX Boo (spectrum 
no. 3; Fig. lA) and /x Cep (spectrum no. 7, Fig. la) are 
included. 

Also included in the figures are ground-based, 
broad-band measurements of the nonvariabie sources 
made by other observers. The agreement is quite good 
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the^oxygen-rich supergiants aOri and ^ Cep. The flux calibration \acs}, other 
relevant data are ®ven in TaWe 2. In this and the other figures, statistical error bars representing ±1 a of the mean arc . shown 
where large enough to plijt. The smooth curve is a 4200 K blaclcbody extrapolated from the 4-8 (im data of Rvjselfet M. and 
**** probable contnbution of the star to the total observed flux from m Cep (spectrum no. 7) descrfibed 
Lnrf open symbols are ground-based broadjband observationa from the following sources: open squares, , Morrison 

Simon, and Dyck 1975; triangle. Low, Rieke, and Armstrong 1973. ‘Ae lower diamond at 
4 Cep spectrum telow it (spec^m no. 7). \b) Tlie observed 16-39 4 m spectra of the oxygen-rich red'igiaut s^ar 
relevant data are »ven in Table 2. The smooth curve is a 2300 K biackbody esttfapolated 
wavelengths and represents our j^timate of the stellar contribution to the spectrum just above. Open square, ground- 
Dasea 20 ftm measurement Qf Momson and Simon 1973. 
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Fto. 1«,— The 16-39 /im spectra of the supeygiant stars VY CMa and IRC + 10420. The flux calibraiion and other relevant data 
are given in Table 2, Broad-band observations by Hagen et al. (diamonds) of VY CMa and by Humphreys et ai, (upside down 
triangles) of IRC + 10420 are included. 
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in all cases except for the 33 urn observation of n Cep 
(Fig. la, spectrum no. 7) reported by Hagen, Simon, 
and Dyck (1975). Their flux level is approximately a 
factor of 2 higher than was observed here, though the 
25 /xm flux levels agree within the errors. A second 
spectrum of p. Cep taken in 1978 January from 16 to 


30 /xm (spectrum no. 8 in Fig. la) agrees quite well in 
shape and flux with the 16-39 /sm si^trum taken in 
1977 June in the region of overlap. Since n Cep is not 
variable in the infrared at 20 fsm (Morrison and Simon 
1973) and has not displayed any large variations in 
visual flux from 1974 to 1978 (Mattel 1978), the 


TABLE 2 


Observino Log and Flux Levels at 20 Microns 


Spectrum, Figure 

Name 

Date 

Notes ‘ 

^20 um 

(W 


(1) 

(2) 

(3) 

(4) 

(3) 

(<S) 


1, la 

.... a Ori 

1976 Nov . 

2, KAO 

1.45 X 10-« 

0.41 

2, 16 

.... RX Boo 

1976 May 

2, KAO 

3.8 X 10-1* 

0.32 

3, 16 

.... RX Boo 

1977 Jim 

2, KAO 

3.9 X 10- 

0.31 

4, 16 

.... RX Boo 

1978 May 

10, KAO 

4.1) X 10- 

0.30 

5, Ic 

X Her 

1977 Jun 

2, KAO 

2.15 X 10-“ 

0.33 

6, Ic 

X Her 

1978 May 

10, KAO 

2.1 X 10-“ 

0.33 

7, la 

.... p Cep 

1977 Jun 

2, KAO 

5.8 X 10-“ 

0.17 

8. la 

.... p Cep 

1978 Jan 

10, KAO 

6.4 X 10-“ 

0.16 

9,1c 

.... PZ Cas 

1976 Nov 

2, KAO 

2.9 X 10- “ 

•CO.l 

10,1c 

.... PZ Cas 

1977 Jun 

2. KAO 

3.7 X 10-“ 

<0.1 

11, Id. 

.... R Cas 

1978 Jan 

10, KAO 

6.8 X 10-“ 

0.22 

12, Id. 

NML Tau 

1978 Jan 

10, KAO 

1.28 X 10-“ 

<0.1 

13, Irf 

OH 26.5+0.6 

1,977 Jun 

2, KAO 

1.08 X 10-“ 

<0.1 

14, Ic 

VY CMa 

1976 Jan 

2, Lear 

6,8 X 10-“ 

<0.1 

15, Ic 

VY CMa 

1976 Nov 

/ 2, Lear 

7.7 10-“ 

<0.1 

16, Ic 

VY CMa 

1977 Sep 

. / 10, Lear 

8.0 X 10-“ 

<0.1 

17, Ic 

IRC +10420 

1976 May 

2, KAO 

2.25 X 10-“ 

<0.1 

18,1c 

IRC +10420 

1978 May 

10, KAO 

2.0 X 10-“ 

<0.1 


* Notes: 2*2 detector 16-39 ;im spectrometer; 10 » 10 detector 16-30 pm spectrometer; KAO m NASA C-141 KAO 91 cm 
telescope; Lear « NASA Lear Jet 30 cm telescope. 

” An entry < 0.1 in column (6) means the star contributes less than 107s of the observed flux at 20 pm. 
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discrepancy is probably not due to source variability, It 
is felt that the measurements reported here arc suffi- 
ciently reliable to rule out the factor of 2 e,xcess for 
/i Cep at 33 jura reported by Hagen et al. This point 
could be settled by an independent measurement of 
the 33 jttm flux level of /. Cep, though this wouid be 
difficult to do from the ground because of the large and 
tirae-variabie opacity of terrestrial water vapor at these 
wavelengths. 

rV. DISCUSSION 

All the stars in this study radiate in excess of what 
would be expected from the star alone in the 16-40 /um 
band (Table 2), The probable source of this emission is 
circumstellar dust heated by stellar radiation. The 
character of this radiation is a smooth, biackbody-like 
continuum at long wavelengths with leveling oflf or 
downturn shortward of ISjum, most evident in the 
cases where the dust emission dominates (e.g,, ju Cep, 
Fig. la; PZ Cas, Fig. Ic; VY CMa and IRC +10420, 
Fig. le; NML Tau, Fig. U). As this downturn is more 
rapid than can result from blackbod>[ emission, it 
indicates that the opacity of the dust grains is not gray 
at these wavelengths. We interpret this as an emission 
feature due to a grain resonance in this region. The 
spectrum of OH 26.5+0,6 (Fig. Id) shows m absorp- 
tion feature at around 18 jum which is due to radiation 
transfer effects in the extremely thick shell surrounding 
the star, as discussed by Forrest et al. (1978). When the 
optical depth is large, a local minimum in flux is 
observed where the grain emissivity peaks due to the 
decrease in dust temperature with distance from the 
star. 

The spectrum of RX Boo (Fig. 16), an M8e semi- 
regular variable, is somewhat different from the other 
stars in that the emission is almost flat out to 20 /xm 
and then decreases. Taking into account the stellar 
contribution at these wavelengths (Fig. 16, smooth 
curve), this indicates the dust emission peaks at 
approximately 19jum as compared with the peak at 
approximately 17.5 jum seen in the other stars. This 
could be due to a different composition of dust grains 
around this star. Interestingly, the 10 /xm spectrum of 
this star (obtained in 1972 [Gillett and Forrest 1978], 
not shown here) is somewhat atypical of M type stars 
in that the excess peaks lon^ard of 9.7 /xm and 
decreases less rapidly beyond the peak. R Leo (Forrest, 
Gillett, and Stein 1975) shows a similar behavior. 

Elitzur (1978) has proposed a mechanism for the 
excitation of the main line maser emission of OH at 
1665 and 1667 MHz which is seen from many cool 
oxygen-rich red-giant stars. In the present study, the 
cool stars with thinner circumstellar shells, such as 
RX Boo, X Her; and R Cas, are typical of this class of 
object. In his model, the ground-state inversion results 
from a nonequilibrium population of upper rotational 
levels due to pumping by a radiation fleld which falls 
more steeply than A « 1/A®. These upper levels are 
populated by photons at wavelengths near A =« 35, 53, 
79, or 119 /xm. Of the stars that have been studied, the 
emergent spectra of the optically thinnest sources with 


the hottest apparent dust temperatures (i.e., with a 
small infrared excess given in Table 1) goes approxi- 
mately as C.C (l/A)* from 30 to 40 /xm. In those stars 
with thick shells the observed flux falls off more 
slowly than this. Thus it is not obvious that radiation 
excitation involving 35 /xm photons as proposed by 
Elitzur (1978) will work. Closer to the star or at longer 
wavelengths the dust spectrum will probably be 
steeper, but this can be offset by a greater contribution 
from the star whose spectrum is probably a blackbody 
in the Rayleigh-Jeans limit and therefore goes at (1/A)\ 

The emission from the dust alone around the stars 
/X Cep, PZ Cas, VY CMa, and IRC +10420 from 8 to 
38 /xm is shown in Figure 2. These stars were chosen 
because of the large contrast between the dust emission 
and the estimated stellar emission at these wavelengths. 

For /X Cep, the stellar spectrum and excess emission 
from 8 to 13.5 /xm was taken from Russell, Soifer, and 
Forrest (1975). Those authors discussed the possibility 
that some of the emission shortward of 8 /xm could be 
from the circumstellar envelope. The stellar spectrum 
in the 16-40 /xm region was assumed to be the con- 
tinuation of the blackbody lit to the S.5^8 /tm data and 
is shown as the smooth curve in Figure la. The 
difference between the observed flux and this black- 
body Is attributed to circumstellar emission and is 
plotted in Figure 2. Integrated over wavelength, the 
circumstellar emission constitutes less than 0.05 of the 
total emission from /x Cep (bolometric excess <0,05; 
Table 1), indicating a relatively thin circumstellar shell. 
The subtraction process is most uncertain near 8 /xm 
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where the stellar continuum can contribute a major 
fraction of the observed Hux. 

PZ Cas was observed from 16 to 39 /xm at the sug- 
gestion of F. C. Giliett (1976), who had obtained an 
8-13 (tm. spectrum in 1975 November showing a very 
large dust emission feature peaking at 9.7 /xm, Using 
the broad-band measurements from the IRC (Neuge- 
bauer and Leighton 1969) avtd AFGL (Price and 
Walker 1976) in conjunction with the 8-13 and 16-39 
’ /xm spectra, it was found that the energy distribution 
could be fitted by a 2500 K blackbody shortward of 
5 nm plus an infrared excess longward of 8 /xm. The 
biackbody was assumed to represent the stellar con- 
tinuum and was subtracted to give the infrared excess 
shown in Figure 2. Integrated over wavelength, the 
infrared excess constituted approximately 0.16 of the 
total flux observed (bolometric excess s 0.16; Table 1), 
indicating a moderately thick circumstellar shell. 
Again, this process is most uncertain near 8 /xm where 
the stellar continuum can supply a major fraction of 
the observed flux. Beyond 10 /xm the stellar continuum 
represents less than 107o of the observed flux. There 1$ 
an additional uncertainty in the ratio of 10 to 20 /xm 
emission because the photometry and spectrometry 
were done at different times and this source appean to 
vary in the infrared (Table 2; Price and Walker 1976). 

For the stars VY CMa and IRC +10420. the optical 
depth of the circumstellar shell is large enough to 
modify the stellar spectrum, and a different procedure 
to estimate the infrared excess was necessary. In these 
f cases the total observed emission was integrated over 
ail wavelengths. In the absence of a circumstellar shell, 
the spectrum of the star was assumed to be a blackbody 
of temperature appropriate for the spectral type 
(6000 K for IRC +10420, 2500 K for VY CMa) with 
the same integrated luminosity as the observed lumi- 
nosity. Any emission above this flux level is certainly 
circumstellar excess. For VY CMa the infrared excess 
integrated over wavelength constitutes at least 0.5 of 
the total flux (bolometric excess ^ 0.5 ; Table 1), while 
for IRC +10420 the bolometric excess was at least 0.88. 
For IRC +10420, the star can contribute at most 2% 
of the flux observed at 8 /xm and beyond, so no sub- 
traction was necessary and the total observed flux is 
plotted in Figure 2. For VY CMa, the star may con- 
tribute as much as 127o to the flux observed at 3 /xm, 
but this fraction decreases to less than 47o at 10 /xm 
and beyond. Since the contribution of the star to the 
observra spectrum is uncertain when the optical depth 
and geometry of the circumstellar shell are uncertain, 
the total observed flux is plotted in Figure 2. It is 
' believed to be primarily due to circumstellar emission, 

but there may be some stellar contribution near 8 /xm. 

From the derived spectra of 8-39 /xm circumstellar 
' emission shown in Figure 2, the similarity of the 

circumstellar dust is apparent. There is a peak of 
emission near 18 /xm with a smooth continuum falling 
toward longer wavelengths, indicating a resonance in 
the grain material near 18 /xm. At shorter wavelengths 
there is the familiar and more prominent 9.7 /xm emis- 
sion feature. The existence of these two resonances is 
; in qualitative agreement with the silicate hypothesis 
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which predicts a feature in the 10 /xm region due to the 
Si^O stretching resonance and a feature in the 20 /xm 
region due to the O-Si-0 bending resonance. 

The smoothness of the spectra indicates the grain 
material is different than well-ordered lunar and 
terrestrial silicates which typically show sharp struc- 
ture, especially in the 20 /xm region (see Knacke and 
Thomson 1973; Zaikowski, Knacke, and Porco 1975). 
The shape of these spectra is quite similar to that 
observed in the Trapezium region of the Orion Nebula 
(Forrest and Soifer 1976; Forrest, Houck, and Reed 
1976), but with a higher apparent temperature. This is 
consistent with the idea that the cool oxygen-rich red- 
giant stars are the source of an important component 
of the interstellar dust seen in Orion. 



V. SHELL MODEL 

a) The Model | 

In order to make a more detailed comparison of the 
observed circumstellar spectra in Figure 2 and the 
theoretical emission, some simple models have been 
considered. In the case where the shell optical depth is 
not large, the theoretical emission is particularl}[ simple 
because each grain is observed without attenuation and 
the radiation field which determines the grain tempera- 
ture falls as 1/R^ beyond a few stellar radii. As there 
are no sharp resonances in the observed spectra, it is 
not possible at this time to unambiguously identify the 
grain material which is responsible for the emission. 
However, the existence of the two features near 10 and I 
20 /xm is consistent with silicate emission, and this will I 
be used as a working hypothesis for the present models. - 1 

It will be shown that using measured grain emissivities 
from a likely silicate material and a realistic density f 
distribution in the circumstellar envelope gives a good 
fit to the observed optically thin spectra. 

Of the laboratory measurements available at present, 
the silicate materials which appear most similar to the 
observed 8-40 /xm spectra are the carbonaceous 
chondrite meteorites Vigarano and Murchison meas- 
ured by Penman (1976) (see Forrest, Houck, and Reed 
1976). This li; physically reasonable because this type 
of meteorite is believed to represent relatively unpro- 
cessed material from the primitive solar nebula. Other 
carbonaceous chondrite materials (Zaikowski and 
Knacke 1975) and an artificially produced amorphous 
silicate material (Day 1974) show a close similarity to 
the 10 /xm feature, but the second peak occurs at 
longer wavelengths than the feature observed here and 
by Forrest, Houck, and Reed (1976). Penman measured i 
the reflectance from polished samples and applied the I 
^mers-Kronig principle to denve the optical con- 
stants from 5 to 40 /im. The emissivities of small I 
spherical grains were calculated using Mie theory. For 
the present comparison, the emissivities of equal parts 
of the Vigamo and Murchison material have been 
taken from Penman (1976); the resultant moss opacity, 
ka, is shown at the top of Figure 3. The peak mass 
opacity of approximately 2800cm‘ g~^ occurs at I 
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Fio. 3.— Model fits to the excess emission from PZ Cas and 
ft Cep (jtoUd llnex) as described in the text. The crosses are for 
a density distribution pt poIR* and the diamonds for a 
density p, - (l/fl’(l - RofRV'^ The assumed inner tem- 
peratures are, from the top, 400 K, 340 K, 1250 K, and 860 K, 
u noted in Tabie 3. The squares are the opacity («a, cm^ g*‘) 
of the meteorite mixture used in calculating the emergent 
spectra. The peak opacity of 2800 cm^ g'^ occurs at 10 pm, 

10 /sm, and a second peak with approximately 1650 
cm^ g~^ opacity occurs near 18.5 pm. 

The flux from a grain of mass mass opacity k^, 
at a temperature Ti, and a distance D away from the 
Earth, will lie 

( 1 ) 

where B/XT) is the Planck blackbody function (W cm"® 
jum"® sr“®). If the circumstellar shell is optically thin 
at infrared wavelengths, each grain will be seen inde- 
pendently, and the total flux observed will be 

( 2 ) 

Thus the emergent spectrum is determined by the 
density of grains in the shell versus distance from the 
star pi(R) and the grain temperature as a function of 
distance from the star T^R). The temperature of a 
grain will be determined by energy balance, 

In the present case will be dominated by therm^ 
radiation, so 

(3) 

Jo 

If the circumstellar shell is optically thin in the near- 


infrared (1-5 ^m), £ia will be dominated by the stellar 
radiation field, which falls as 1/F® more than a few 
stellar radii from the photosphere, Therefore £,a * 
£o{Ro/RT> where is the power input at a distance 
R m Bq from the star, 

These equations thus define the run of temperature 
versus distance from the star. In order to evaluate 
equation (3), it is necessary to know the grain opacity 
fts at ail wavelengths. From 5 to 40 pm, the opacities 
from Penman (1976) shown in Figure 3 were used. 
Shortward of 5 pm, it was assumed the optical con- 
stants were constant versus wavelength and the par- 
ticles were small so 

/ij Ml 

» 

Beyond 40 Mm, it was assumed that the imaginary part 
of the dielectric constant fell as 1/A so that 

1 40 pm V 

■* — j ‘ 

At the temperatures of most interest in the circum- 
steilar shells between 100 and 1000 K, more than half 
of the Planck emission falls between 5 and 40 ptm, so 
it is believed these approximations are not critical to 
the results. Numerical integration of equation (3) 
showed that Fom oc T* * from T ■ 160 to 1600 K and 
Boat « below 160 K. 

The density distribution in the circumstellar shell is 
less certain. The mass-loss rate from a star is given by 

Air - 4rrR^p(R)v(R) , (4) 

where p{R) is the total density and v(R) is the outflow 
velocity. Since both the gas and the dust mass-loss 
rates have been observed to be fairly constant from 
stars of this type (see Forrest, Gillett, and Stein 1975; 
Weymann 1963) and the circumstellar absorption lines 
(Deutsch 1960; Weymann 1963) and OH emission 
lines (Wilson and Baitei; 1972) indicate the outflow 
velocity is constant with time, it is e,tpected that 
p oc 1/R4 in the outer parts of the envelope. Since dust 
grains are not expected to condense where the density 
is very low and the dust is momentum coupled to the 
gas (Gilman 1972), one also expects the dust density 
Pi oc 1/R® in the outer parts of the envelope. Close to 
the point where dust grains begin to condense, neither 
the velocity nor the fractional amount of dust is well 
known. Fortunately, it will he shown that at the wave- 
lengths of interest beyond 8 Mm, the shape of the 
emergent spectrum does not depend strongly on the 
density distribution near the region of condensation. 
The point at which dust begins to condense will be 
denoted Ro ^ r«, where r« is the stellar radius,. Theory 
(Jones and Merrill 1976) and observation (TAppala ei 
al. 1974) indicate Ro/r^ ^ 3 for the type of stars con- 
sidered here. For the purposes of this calculation, it 
has been assumed that all the dust condensation 
occurs at Ro (see Menietti and Fix 1978). Then if 
radation pressure on the grains determines the flow 
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velocity, the dust velocity wili be of the form 

Vd(Ji) « vM + Vn{i - , (5) 

where Vd(Mo) is the initial velocity and Vi(Rn) + v„ is 
the final outflow velocity, Two cases have been con* 
sidered: v^iRo) » which leads to a pi oc IjR? shell 
structure, and " 0, for which 


■> R4Ry>^ 

IS infinite at Ro but quickly approaches 1/i?* beyond 
2 R(). 

In a preliminary run of the model with Pi <x l/i?“ 
and Ttt » 500 K, an interesting effect was discovered 
when the model v/as run out to outer limits of 10 i?o» 
100 Ro, and 1000 Rq, Though 997o of the total flux 
comes from within 100 Ro, this region contributes 
only 907e of the flux at 40 fitn ann approximately half 
the flux at 80 urn, Even though the temperature is low 
(~ 35-76 K) in the outer layers of this shell, the mass 
is large, and this region dominates the emission at 
longer wavelengths. Farther from the star than 
approximately 1000 Rq, two effects are important in 
determining the emergent spectrum beyond 100 /xm. 
First, heating by the interstellar radiation field will set 
a minimum temperature of the grains. With the grain 
emissivity used here, a grain radius of O.l /xm, g»b. * I, 
and the interstellar radiation field given by Allen 
(1973), the minimum temperature is approximately 
30 K in the galactic plane, f’econd, when a mass of 
interstellar matter equal to that in the circumstellar 
shell has. been swept up, the structure of the shell will 
be significantly different than the /> oc l/R!^ assumed 
here. For a space density of I H atom cm”®, an out- 
flow velocity of lOkms”^, and a mass-loss rate of 
10"® yr"S this happens after approximately 3 x 
10® years at a distance R > 10* /f, from the star. 
Since we are interested in the emergent spectrum from 
5 to 40 ixm, these effects will not be important and it is 
sufficient to integrate out to 1000 Rq. 

With the density in the shell fixed, the only free 
parameter in the model is To, the dust temperature at 
the distance Ro from the star where dust forms. This 
temperature has been varied, in order to fit the observed 
ratio of Fio umIRia m fof the shells around PZ Cas and 
jj. Cep (Fig. 2). The temperatures are given in Table 3 
and the resulting spectra are compared to the observed 
spectra in Figure 3. Also given in Table 3 is the tern- 
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perature a thin isothermal shell would have in order 
I’o give the observed ratio 

b) Discussion 

From Figure 3 it is seen that the models give a good 
fit to the observed spectra. In particular, the shape and 
position of the 10 and 18 /xm emission features are 
reproduced, This provides further support to the 
hypothesis (Woolf and Ney 1969) that silicate grains 
are responsible for much of the excess emission from 
oxygen-rich red-giant stars. It also indicates that the 
grains around these stars are probably made of a dis- 
ordered, amorphous form (if silicates similar to that 
found in the carbonaceous chendrite meteorites. 

From Figure 3 it is seen that there is little difference 
between the emergent spectra for the two density 
distributions considered here. This shows that the 
emergent spectrum is not particularly sensitive to the 
density distribution near R Ro, where it is least well 
known. Thus it is feit that the discrepancies between 
the model and observed fluxes shortward of 9 and 
longward of 20 /xm are not due to density effects alone. 

Short of 9 /xm, the models predict more flux than is 
seen. This could be due to either lower emissivity of 
the circumstellar grains at these wavelengths or an 
overestimate of the stellar continuum which was sub- 
tracted to derive the excess emission. For both PZ Cas 
and /X Cep, the 8-10 /xm spectrum had the approximate 
shape of the model spectrum before subtracting the 
stellar contribution, This indicates that the circum- 
steilar silicates probably have a sharper feature than 
the silicates measured by Penman (1976). However, 
some of the observed 6-8 /xm continuum could be due 
to the hot grains in the inner part of the shell. A 
measurement of the size of the emitting region at 
different wavelengths could clarify this point. 

Longward of 20 /xm, the models predict less flux 
than is observed. This indicates that the emissivity of 
the circumstellar grains falls less rapidly than the 
silicates measured by Penman (1976). The observed 
spectra would be fitted by an emissivity which fell as 
l/A® from 20 to 40 /xm rather than the l/A®*® depend- 
ence of the Vigarano and Mirchison mixture used here. 

The values found here for the dust emissivity and 
dust temperatures in the shell are consistent with the 
model or Elitzur (1978) for producing main line OH 
maser emission mentioned earlier. However, a detailed 
comparison for specific stars is necessary to sec 
svhsther the shell optical depths are large enough to 
produce a sufficiently intense incident flux with a steep 
enough spectrum in the masing region. 
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TABLE 3 

Model Dust TEMPERATxniES 



Fio umIFii urn 

(observed) 


To 


Star 

Thin Shell 

Pi - PJR’‘ 

Pi » (Po/i?®(l - 

PZ Cas 

2.95 

270 K 

400K 

340 K 

M Cep 

6.9 

470 K 

1250 K 

860 K 


620 
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VI. CONCLUSIONS 

The study of the 16-39 fitn emission from oxygen- 
rich red-giant stars and IRC +10420 with excess 
emission at lOjum has revealed: 

1) The sUrs show excess emission in the 16-39 litn 
region in the form of a broad hump peaking near 
IS (im and falling smoothly to longer wavelengths. 

2) Except for possibly RX Boo, the emission from 
these stars is quite similar in character and qualitativeiy 
like the emission seen in the Trapezium region of the 
Orion Nebuia (Forrest, Houck, and Reed 1976). This 
indicates the grain material are quite similar in these 
objects. 

3) Lack of thither structure in the 20 nm region 
prevents an unambiguous identification of the circum- 
steiUu material. The observed spectra are consistent 
with emission from small grains of silicate material 
similar to the carbonaceous chondrite meteorites in a 


diflfuse circumstellar envelope. Comparison of the 
spectra with a simple model indicates that the grain 
emissivity falls as 1/A® from 20 to 40 nm and suggests 
that some of the observed 6-8 urn continuum could be 
due to emission from hot grains in the inner part of the 
circumstellar envelope. 

4) The spectra of fi Cep do not show the 33 fxm 
emission feature reported by Hagen el al. 
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THE GALACTIC CENTER! 16-30 MICRON OBSERVATIONS AND THE 

18 MICRON EXTINCTION 
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ABSTRACT 

Low resolution 16-30 /im spectra of three regions In Sgr A West have been obtained, and the 
area has been scanned in two colors, \9 urn and 28 /sm, with 30" resolution. The J8 urn silicate 
feature is clearly seen in absorption as a broad depression present in all three spectra. Two- 
component models of the continuum have been constructed and fitted to the spectra. The most 
plausible models assume silicate grains for both components and derive a 19 /an optical depth for 
^e absorbing dust of —2. Comparison with far-infrared data suggest that the absorbing component 
is in fact made up of two parts: one local to the galactic center, and the other distributed along the 
10 kpc line of sight, with the local component constituting roughly one-third of the total. Extension 
of the models to \< 15 /tm shows that additioiud sources of emission are necessary to explain the 
observed flux at shorter wavelengths, and that the ratio of optical depths in the 9.7 and 18 /tm 
siUcite feature niay be as large as 2=3. The total dust mass column density in the line If sight is 
estimated to be —1.4 X 10"^ g cm**^, The local component is within —2 pc of the galactic center 
and has a mass of —20 Mq, 

No 18.7 /tm [S lu] emission has been detected, which Is consistent with normal S m/H n 
abundance and an 18.71 /tm dust extinction optical depth of 2.1 or more, in agreement with the 
most plausible continuum models. Indirect arguments suggest that the S/H abundance is no more 
than a factor of 2 higher than in the solar neighborhood. G^mparison with the Brackett-y flux 
suggests that the S/H abundance is lower than this. 

Subject headings: pAaxieat Milky Way — galaxies: nuclei — infrared: spectra 


I. INTRODUenON 

Within the galactic center, ionized gas, stars, and 
heated dust are in close proximity to each other. The 
ionized gas is observed in the thermal radio continuum 
source caUed Sgr A West, in many infrared and radio 
recombmation lines of hydrogen and helium, and in 
infrared fine-structure lines such as [Ne n] 12,8 /tm and 
[0 m] 88 /im (Eken et al. 1975; Neugebauer et aL 
1978; Mezger and Smith 1976; Lacy et al. 1979; Dain 
et al. 1978). Sgr A West is also a strong iiifrared 
continuum source. In the near-infrared (X;g3/tm), the 
flux is mainly due to the late-type stan cluster^ at the 
core, attenuated by the intervening interstellar dust 
(Becklin and Neugebauer 1968). From 3 /tm to 1 mm, 
the strong infrared continuum is produced primarily by 
thermal radiadon from clouds of dust grains. In the 
mid-infrared (3 /tm'<\:S20 /tm), ground-based high- 
resoludon maps show several compact continuum 
sources (BeckUn et oL 1978a; Rieke, Telesco, and 
Harper 1978). Some, for example IRS 7, are identified 
at late-type supergiaut stars, but others, such as IRS 1, 
are not well understood. Some evidence exists that they 


are compact clumps of ionized gas with associated dust 
(Lacy et al. 1979). The 8-13 /tm spectra of ail the 
sources m Sgr A West are dominated by deep 9.7 /tm 
absorption features which are attributed to attenuadon 
by cold foreground silicate grains of the background 
continuum from warm dust (Woolf 1973; Willner 197Q. 
At A— SO /tm, available maps have lower spatial resolu- 
don, and only two sources are identified (Harvey, 
Campbell, and Hoffmann 1976; Rieke, Telesco, and 
Harper 1973). One, Sgr A West (NE), is comcident 
with the thermal radio source and the compact mid- 
infrared iiources. The other, Sgr A West (SW), has a 
mt«ch colder spectrum. Although NE and SW are ap- 
proximately the same brightness at 50 /tm, at 10 /tm 
SW is much dimmer, and to date only a single compact 
source has been detected in it (Lebofsky 1979). 

In this paper we will address several quesdons of 
current interest First, what is the amount and distribu- 
don of the dust both within the galacdc center and 
along the line of sight? What is the physical tempera- 
ture of the grains? What can we conclude about the 
emissivity spectrum of the grains, and what does this 
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Imply about the eompositfon of the grains? As for the 
Ionized gas, what is its ionization structure and what 
are the elemental abundances in It? What is the nature 
of the infrared continuum sources? 

To attack these questio;ts, we have obtained low* 
resolution 16-30 pm spectra of three positions in the 
Sgr A West region. In addition, we have mapped the 
area in two colors (19 and 28 pm) with a spatial 
resolution of 30" (1.5 pc at 10 kpc); these results are 
presented in | II. In § III, upper lindts to the 18.71 pm 
(S 111 ] flux are given and interpreted in terms of the dust 
extinction and the ionization structure and elemental 
abundances of the ionized gas. Section IV discusses the 
continuum emission and proposes models to explain it. 
An 18 pm absorption feature Is observed and identified 
as a resonance of silicate minerals present in dust 
grains. Section V examines the results of the continuum 
models, determines the consequences of extending the 
models to other wavelengths, and discusses the relative 
strengths of the 9.7 pm and 18 pm silicate features and 
the shape of the silicate spectrum for A > 16 pm. In § 
VI, we present our conclusions. 

II, OnSBRVATIONS AND RESULTS 

Observations of the galactic center were made during 
two flight series on NASA aircraft flying at altitudes of 
more than 12 km. In both series the instrument used 
was a iO*channel liquid helium cooled grating spec* 
trometer described by McCarthy, Forrest, and Houck 
(1979). The first observations were performed in 1977 
September on the Lear Jet, and the second set in 1978 
May using the Kuiper Airborne Observatory (KAO), 
Relevav'xt observational pora<neters are summarized in 
Table 1. Standard beam-switching and data-reduction 
techniques were used. In the calibration spe.:?Aim of 
Mars from the 1978 May series, a narrow water vapor 
feature was noted at 29.3 pm and was suppressed with 
a 15% correction to the calibration. 

A continuum spectrum from 16 to 30 pm was ob- 
tained in the Lear Jet series. The 217 ffeld of view ("big 
beam”) was large enough to contain the entire Sgr A 
West source, which has a FWHM diameter of —30" in 
the 16-30 pm range. During the 1978 May KAO 
series, 16-30 pm spectra were obtained of both far- 
infrared sources, NE and SW, at a resolution of 30". 
For the SW spectrum, the beam was nearer the position 
of Harvey, Campbell, and Hoffman (1976) than that of 
Rieke, Telesco, and Harper (1978), but at 24 pm this 
source was not clearly separated from the peaic (NE). 
The reduced spectra are shown in Figure I. In both of 
the KAO spectra there are systematic enors at \ > 20 
pm. The fint, type, so-called scalloping, appears as a 
wavelike variation with a period of 2 pm and a peak-to- 
peak amplitude of 10- 15% of the flux. This appears to 
be due to phase differences between channels, aggra- 
vated by the high chopper frequency (35 Hz). Also, two 


TABLE I 

OlSeXVAtlOKAL Paramstsm 


Fuoiir Semes 

PAXAMmx 1977 s«pt 1978 M*y 

AlfCfflt/t IfCAT 

C-Nl (KAO) 

Telescope* 30 cm 

91 cm 

B<«ndiimctcr(FWHM).,. 2'.7 

30' 

Chopp«Tlbrow., 8' 

90* 

Spectral r**olutton (FWHM) 0.5 n m 

02ftm 

SUndwd spectrum. .Moon 

Mm 

Fluxsundard,,,,,,,.,,.,. Mars 

M*n,aHef,IRC+10420 


peaks at 26.5 and 28.2 pm show the effect of Ino' p- 
pletely canceled water vapor features. None of these 
features appear to be real properties of the source 
spectra. 

Several comments may be made here: First, all the 
spectra have a broad depression centered at —19 pm. 
We attribute this feature to absorption of radiation by 
a foreground cloud of cold silicate grains of the back- 
ground radiation from hot dust. Below we report on 
more specific models of this continuum. Second, a 
comparison of the Sgr A West (NE) and Sgr A West 
(2!7) big beam spectra shows that the central 30" con- 
tains from one-half to one-fourth of the total flux from 
the entire galactic center region, indicating that the 
source is not much larger than the 30" beam. Further- 
more, since the color temperature of the big-beam 
spectrum is lower than that of the NE source, the outer 
regions of the source must be colder than the center, 
and the source diameter must increase with wave- 
length. Both of these effects are confirmed by our data. 
Ground-based 21 pm maps (Becklin et al. 1978a) show 
that the galactic center contains a number of compact 
sources in addition to an extended component, but 
they contribute only about 20% of the flux in a 30" 
beam. Therefore, these compact sources are not vital to 
our understanding of the present spectra. Third, the 
SW source and the big beam spectra have virtually the 
same shape, suggesting that the Lear Jet spectrum 
might be modeled as a source similar to the SW one 
but filling roughly one-tenth of the big beam’s solid 
angie. Fourth, none of thi); spectra show evidence of [S 
in] 18.71 pm emission from the ionized gas, althou^ 
this is a common characteristic of other galactic H u 
regions (McCarthy, Forrest, and Houck 1979; 
McCarthy \^80). 

In 1978 May we also observed 14 positions at two 
different wavelengths (18.9 and 27.8 pm) and con- 
structed a map of the Sgr A West region. The points 
are arranged in a rough cross with the long arm ap- 
proximately parallel to the galactic equator and with 
the NE source located at the intersection of the arms. 
Individual points are separated by —12", which is 
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Fio. I. — Reduced spectra of the galactic center. The top spectrum was taken with a 2.7' beam and includes the entire source. The 
middle spectrum wu centered on the peak 24/im source in a 30* beam, which approximately coincides with the NE source. The bottom 
spectrum, also with a 30* beam, is of a point near the SW source. The curves are model fits to the data. The model consisu of optically 
tUn dust obKured by a cold dust cloud with optical depth of order unity. Both clouds are assumed to consist of silicate dust grains. 


slightly less than a beam radius. The chopper throw 
was sufficiently large that the flux in the reference 
beams was ne^gibly small. In the coune of the map- 
ping we measured the NE source three times (positions 
1, 4, and 14), and these results can be used to gauge the 
effects of guiding jitter. At the position of the NE 
source the graphs show an average of these three mea- 
surements. 

Beam positioning on the KAO was done with a p-^ 
offset stage, and by peaking up on the 24 /im con- 
tinuum source. The offset guide star was an anony- 
mous star which is located 135" from the NE source 
(1950 position 17'*42“20.*3, -28"58'41" from Palomar 
plates). Inaccuracies in the computed field rotation 
angle dominate the relative positional )incertainties 
which are estimated at —3". The uncertainty in the 
absolute positions depends on the position of the pMo 
star (known to ~± 10") and the absolute error , , 
rotation angle. G>mputing the rotation angle has 
a constant source of trouble in the past, and it is 
difficult to estimate its uncertainty. Since the peak flux 
of the 19 nm map agrees well wi^ ground-based mea- 
surements, we have corrected the derived position of 
the NE source to agree with IRS 1. This correction was 
then made to ail the other map positions and to the 
positions of the full spectra. It corresponds to an error 
in the calculated field rotation angle of 5% or '«12" on 
the sky. It should be noted that the beam position of 


the SW spectrum does not agree well with the positions 
of Harvey, Campbell, and Hoffman (1976) or Rieke, 
Telesco, and Harper (1978). At best we were pointed at 
the edge of the source and included about half of its 
solid angle. However, we did include most of the weak 
and extended 10 jam emission associated with this 
source. 

The maps are shown in Figure 2 and the raw data 
are listed in Table 2. The fluxes have an effective 
bandwidth of 2 pm since they were derived from an 
analog average of all 10 detector signals. Again, the 
quoted flux uncertainties are one standard deviation of 
the mean as derived from the scatter of several integra- 
tions. 

In Table 2B and Figure 2, two quantities are derived 
from the raw fluxes and graphed. The first is the color 
temperature, and the second is the ratio of the observed 
flux to the flux expected from a blackbody which has 
the observed color temperature and which fills the 
beam. This second parameter is a very crude measure 
of the source’s apparent optical depth. One must be 
careful about simplistic interpretations of these quanti- 
ties, however. Since there is an obvious spectral feature, 
the spectrum is not that of a blackbody, and the color 
temperature does not correspond to the temperature of 
dust in the source. Furthermore, it is often the case in 
H U regions that a distribution of temperatures is 
necessary to explain the spectrum adequately,.., 



SBL * 

• • 

fze» 


Sgr A West 


18.9/i.m Continuum 


29.8 /im Continuum 


- 28 * 59 ' 


St 


i9yu.m/28fim 
Color Temperature 


17"42"*30* 


■29*00'H«) 


I X 10^ 

Bx(Tc)iiB 

^ . . I ■ 

17^42" 30* 


Right Ascension (1950) 


Fio. Z— (a) FeatUKS is Sgr A West Orel*, IRSl 20 ftm peak (Becklin «f aL 1978a). Soiid lint, Outensott contour of 20 /tm map of 
Bccklin tt oL (1978a). Triangle, sguare, and dashtd line, 53 peolu of NE and SW souicct, and 53 fim half nuuc contour (Harvey et aL 
1976). Dotted line. Half max 5 GHz contour of WORST map (Been «r oL 1975). (b) The 18.9 pm continuum of Sgr A West The units are 
10" '* W cm~^ um~ and in these units a typical error is ±0.3. Beam size— 30*, effective bandpau— ZO pm. (c) The 29.8 pm continuum. 
Units are 10~ '* W cm~^ pm~ ' , and a typi^ uncertainty is ±02. (<f ) The 19 pm/28 pm colw temperature map of Sgr A West derived 
from Figs. 2b and 2c. Units are keivins, and typical erron are ±3 K. (<) A map of the ratio of the observed flux to the flux from a 
biackbo^ of the derived color temperature which fills the beam. This quantity is a crude measure of the apparent optical depth and is 
dimensiocless. The plotted numben are to be multiplied by 10~^, and a typical eiror is ±5%. 









TABLE 2 

SAorrrAJUus A West Map 


A. Observed Quakitties 


PosmoN 

No. 

PosmoN 

Continuum Flux 
( 10-1* Wcm”^/im"‘) 

Notes 

R.A. (1930) 

Decl. (1950) 

at 18,9 ;im 

at 27.8 nm 

1 

17‘42"30?1 

-28“5920' 

16,1*0.2 

14.9*02 

• 

2 

31.3 

38 36 

4P*0.1 

6.1*02 


3 

30.8 

3908 

10.9±0.6 

11.7*02 


4 

30.2 

39 18 

15.8±0.2 

14.9*0.2 


5 

29.6 

3928 

12.6*0,1 

10.7*0.4 


6 

28.9 

5936 

5,3 *0.8 

7.3*0.1 


7... 

28.1 

39>43 

2.6±0.2 

32*0.1 


8 

27.4 

5949 

0.9*0.1 

3. 4 *0.1 

b 

9 

26.6 

39 53 

0.7*0.1 

2.1S0.1 


10 

31.7 

3844 

UiO.2 

32*0.4 


11 

30.2 

5918 


• * • 

itd 

12 

28.7 

39 12 

9.1 *0.3 

9,3*0.1 


13 

29.4 

5913 

14.9*0.3 

14.1*02 


14 

30.2 

39 18 

13J*0.3 

14.7*0,1 


15 

30.9 

5921 

S.8*0.2 

92*02 


16 

31.7 

59 24 

4.1*0,! 

4.6*02 


17 

27.8 

3909 

4J*0.2 

32*02 


NE avg, , . i 

30fJ 

3920 

15.?±0,4 

14.8*0.1 

■tC 

(NE) 

30*2 

39 16 

14.5*1.0 

14.3*1.0 


(SW) 

28“3 

39 39 

3J*0.5 

3.6*0.3 

b.« 

2'.7 

» « • 

... 

34*2 

31*4 

• 

B. Derived QuAminEs 

■Position 

19 Mtn/28 fim Color 

F, 



N^. 

Temperature (K) 

X 10“ 

Notes 

1 

133 ±1 


6.0*0.! 

ft 

2 

113±2 


4.6*0.1 


3 

123*3 


6.3 *0.3 


4 

131*2 


6.4±0.1 

ft 

5 

140*3 


3 .3 ±0.2 


6 

109*8 


7.2*0.8 


7 

94±3 


11*0.4 


8 

74*3 


31*0.9 

b 

9 

81*3 


11*0.6 


10 

88*6 


9.7*0.S 


11 

* « • 


,, 

ite 

ft,d 

12 

123*2 


5.0*02 


13 

131*2 


6.1*02 


14 

130*2 


6J*0.1 

ft 

13 

123*2 


4.8±0.1 


16 

120*4 


2.9*02 


17.... 

116±3 


3.8*02 


NE avg. . . . 

131*2 


6.4*02 

ft.C 

(NE) 

127 


7.1 


ft*« 

(SW) 

99 


.8,9 


b,ft 

2'.7 

106 


1.3 


0 


Notes to Table 2 : 

•SgrA West(NE). 

‘’Near Sgr A West (SW). 

“Average of # 1, ?r4, and # 14. 

‘^No fflcasurements made. 

“Complete spectra. 

is the observed flux, is the Planck function for the' derived color temperature 7^, and 
Qj is the solid angle of the bea^ 
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The scans through the source show that, as suggested 
by the full spectra, the color temperature drops as the 
beam is moved off the NE source, but the ratio of the 
flux to the equivalent blackbody flux is nearly constant 
except for an increase near the SW source. This il- 
lustrates that the SW source has a cold spectrum. 

Our spatial resolution is barely good enough to re- 
solve Sgr A. On the NE-SW arc the apparent FWHM 
at 19 lita is ~40" while it is 48" at 2$ /im. Similarly, on 
the E-W line the source shows a FWHM of roughly 
36" at 19 /xm and 40" at 28 /xm. Along both axes the 
source is wider at 28 /xm than at 19 /im, and it is 
slightly elongated along the galactic plane, in agree- 
ment with other maps. It is also apparent that the 
position of peak flux lies between positions 1 and 13, 
slightly to the west of the NE spectrum position. 

ni. [S ra] UPPER LIMITS 

No unresolved emission lines are .<cen in any of the 
three spectra. This is surprising, because at this resolu- 
tion the 18.71 /xm [S lu] line is commonly observed in 
galactic H n regions (McCarthy, Forrest, and Houck 
1979; McCarthy 1980). In Table 3 we present 3 a upper 
limits to the line flux. Using a simple method described 
by McCarthy, Forrest, and Houck (1979), we have 
cdculated the expected line fluxes. This model uses the 
observed thermal radio continuum flux as a measure of 
the mass of ionized gas and assumes that the S m/H n 
abundance ratio observed, in Orion (Peimbert and 
Torres-Peimbert 1977) may be applied to Sgr A West 
Furthermore, emissivities are calculated in the low den- 
sity limit (/I, <6000 cm”^), and 18.71 /xm extinction is 
neglected. As shown in Table 3, the predicted line 
fluxes are 3 to 8 times higher than observed. 

There is evidence for clumps of n,-~10^ cm“^ em- 
bedded in the lower density (n,~-10* cm"^) ionized 
gas (Lacy et al. 1979; Ekers et cd. 1975; Rodriguez and 
Chaisson 1978). However, these clumps contribute only 
about one-fourth, of the total radio flux, and they 
cannot explain the low observed [S m] fluxes. 

The 18.71 /xm extinction to the ionized gas b likely 
to be considerable. Continuum models, presented be- 
low, give t( 19 /xm) up to 2.4, which is high enough to 
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bring the predicted fluxes into agreement with the 
observed upper limits. Therefore, our data are con- 
sistent with normal S iii/H u abundance and large dust 
extinction [t( 18.7 1 /xm)s2.1). 

Compared to other galactic H ii regions, the ioniza- 
tion structure of Sgr A West is peculiar, and the abun- 
dances of heavy elements may be higher than near the 
Sun, Within large uncertainties, the He ii/H ii ratio is 
normal (0.095 ± 0.040), but the line profiles suggest that 
He/H might be abnormally high (Mezger and Smith 
1976). From observations of [Ar ni] 8.99 /xm, [Ar u] 
6.98 /xm, [0 ui] 88 /xm, and [Ne n] 12.8 /xm, there is 
clear evidence for an unusually cool radiation field 
(Lacy et al, 1979; Willner et al. 1979; Dain et al. 1978), 
Furthermore, the [Ar n] and (Ne u] fluxes suggest 
overabundances of these elements. Under these condi- 
tions of partial helium ionization and low Ar m/Ar ll 
ratio, theoretical models suggest that S fn/S^0.5 
(Balick and Sneden 1976; Pipher, private communica- 
tion). Since normal H n regions have S lu/S near unity 
(McCarthy 1980), the S/H abundance ratio can be no 
more than a factor of 2 above the value in the Sun’s 
neighborhood, as long as these theoretical H il region 
models are correct. 

The above estiimate of S m abundance in the galactic 
center depends ott an estimate of the amount of ionized 
gas withm our 30" beam based on the radio map of 
Ekers et al, (1975). The data of Neugebauer et al. 
(1978) on the eniission of Brackett-y (By) from the 
ionized gas indicate that the gas may be more cen- 
trally condensed than indicated by the radio map. 
Using the By results of Neugebauer et al, (1978) in a 
32" beam to c-stimate the amount of ionized gas in our 
beam results in an estimate n (S m)/n (H)^ 1.5 X 10“* 
e}ip(T(18.7 /xm)— t( 2.16 /xm)} for an electron tempera- 
tijire r,^5000 K. As discussed by Neugebauer et al. 
(1978), the extinction to the ioniz^ gas in the galactic 
center is around T 2 .u ^b>^ 2.5. There may be additional 
Ic>cal extinction. Then, if we estimate rnj 
’’’ 9.7 the abundance is n (S ni)/n 

(H)^0.4X10■^ This is less than one-quarter me 
“cosmic” abundance of sulfur (Allen 1973), and indi- 
cates either a low sulfur abundance or a dominance of 


TABLES 

[S m] Upfer Loan 


SncniuM 

Continuum 

(10-'*Wcm->/xm-‘) 

.^(SGHz) 

(Jy) 

Line Flux 
(I0-“Wcm-*) ' 
Observed Predicted 

Observed 

Predicted 

Sgr A West 

3.40±0.04 

26 

<1.9 

5.8 

<0J3 

Sgr A West (NE) . . . 

1.475 ±0.006 

6 

<0.16 

133 

<0.12 

Sgr A West (SW)... 

0.3a7±0.004 

4 

<0.22 

0.89 

<0.25 


Note. — S GHz radio contmuiim fluxea are based on the WORST map of Ekers et al. (1975). Observed line 
fluxes are 3« upper limits. 
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the S II ion. This is surprising in view of the high 
abundance of Ar ii and Ne ii which have similar 
ionization potentials. A low sulfur abundance seems to 
be indicated. Because the extinction at 4.0S /tm is 
closer to that at 18.7 fim, a better estimate of the S in 
abundance could be made if Ba measurements were 
made in a 30" beam. 

IV. CONTINUUM E.MISSION AND .MODELS 

The most striking feature of the continuum spectra is 
the 18 /im depression, which we attribute to absorption 
by sold silicate grains along the line of sight. The 18 
nm silicate feature is commonly observed in oxygen-rich 
stars (Forrest, McCarthy, and Houck 1979) and in H n 
region/moleculor cloud complexes. It is due to O-Si-O 
bending motions in the grains. In both types of sources, 
the feature is seen over a range of optical depths 
running from r<l in emission to 1 in absorption. 
From models of the circumstellar shell emission of 
oxygen-rich stars an empirical emissivity function 
for interstellar silicate grains has been found (Forrest, 
McCarthy, and Houck 1979). This function is plotted 
in Figure 3 and discussed in more detail below. 

To explain the 18 /im feature as absorption, we need 
to assume at least two components to the dust: one to 
do the emitting and the other the absorbmg. The emis- 
sion is from a hot, optically thin cloud of dust grains at 
the galactic center. The absorption is by cold silicate 
grains with 18 /im optical depth of order unity between 
us and the galactic center. These grains may be inter- 
stellar or they may be local to the galactic center or, 
most likely, a little of each. Using simple radiative 
transfer theory, the flux received at the earth is: 

F-B(r,) [1 -cxp(~T,)] exp(-T 2 ) a 
+ B(7’2)[l-.e,xp(-T2)]a. 



Fio. 3.— Normalized einiuivity assumed for the Intcistciiar 
silicate grains. Tlus emissivity fits the spectra of circumstellar 
shells. Beyond 19 ^m, it follows a law. 


Here Tj and Tj arc the temperatures of the hot and 
cold grains, respectively; tj and tj are the wavelength- 
dependent optical depths; 11 is the solid angle of the 
beam; and 3(T) is the Planck blackbody function. Our 
model explicitly includes the emission from the cold, 
absorbing component since far-iiifrared observations 
show that the cold dust can be responsible for a signifi- 
cant fraction of the emission at 30 fim. Note that this 
model also assumes that both components fill the beam, 
an assumption which should be more reliable for the 
KAO spectra. 

This is the simplest model which can produce the 
observed absorption feature, and it is easy to criticize 
it. Due to the observed nonuniform distribution of 
late-type stars and ionized gas, it seems clear that there 
will be a continuous distribution of dust temperatures. 
The split mto two components is then only an ap- 
proximation which will work to the extent that two 
temperatures dominate over all the others. We shall 
see, however, that the model does remarkably well in 
explaining the 16-30 /im data, and with only minor 
modification it can be extended to longer wavelengths. 

Tht» emission of the hot dust is not known, but we 
can make reasonable guesses which will span the possi- 
bilities. Three alternatives have been considered. In the 
first, the grain emissivity was taken to be constant, 
producing a blackbody spectrum for the hot dust cloud. 
This is rather implausible on physical grounds, since 
either the grains would have to be large (radius >5 pm) 
compared to typical interstellar dust or the cloud itself 
would need to be dense enough to be optically thick. 
The second alternative used the previously mentioned 
empirical silicate emissivity function for the hot dust, 
resulting in an emission spectrum similar to that of the 
Trapezium. For \> 19 pm the emissivity falls as 
and for \<19 pm it follows the emissivity of the 
carbonaceous chondrite meteorites Vigamo and 
Murchison as measured by Penman (1976). This curve 
was found to be a good fit to the silicate emission 
spectra of circumstellar shells (Forrest, McCarthy, and 
Houck 1979). Since the temperature of the circumstel- 
lar dust is not known independently, the shell spectra 
do not uniquely determine t(9.7 pm)/T(19 pm). The 
present models are less sensitive to temperature effects 
and may allow us to determine t(9.7 pm)/r(19 pm). 
For a fi^ alternative, we considered a grain emissivity 
which was proportional to This produce a smooth 
spectrum qualitatively si milar to a blackbody spectrum 
of higher temperature. One would expect such an emis- 
sivity from sn^ grains (radius<5 pm) with a constant 
complex index of refraction. Such a grain emissivity is 
often mentioned in the literature, but it is not yet well 
justified observationally. On observational grounds, the 
silicate emissivity must be considered the most likely, 
but the others should not be ignored. 

The two-component model has four free parameters, 
r„ri,r:, T;, and the best fit Values were found by 
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TA2LE4 


CoKHNUVM Fits 


SpecnwjM 



Hot Dust 



Cold Dust 


Type*’ 

T{K) 

T(19^m) 

Type" 

n/C) 

T(19pm) 

SgrAW(NE).. 

0.68 

Si 

184 

9.0X10-’ 

Si 

50,4 

2.1 

S|rAW(SW).. 

1.59 

SI 

184 

2.SX10-’ 

Si 

49.9 

2.4 

SirA(2'.7) 

0.78 

Si 

154 

1.6X10-’ 

Si 

41.8 

2.4 

S|rAW(NE).. 

0.68 

1/X 

149 

6.8X10“’ 

Si 

53.8 

0.84 

S|rAW(SW).. 

1.39 

lA 

139 

1.7X10-’ 

Si 

54.9 

0.78 

SgrA(2'.7) 

0.78 

lA 

125 

1.3 X 10“’ 

Si 

42.8 

1.09 

SgrAW(NE)., 

0.39 

BB 

187 

2.6X10-’ 

Si 

49.3 

0.9 

S*rAW(SW).. 

1.48 

BB 

189 

6.0Xl0-‘ 

Si 

52.4 

1.1 

S|rA(2'.7) 

0.76 

BB 

150 

5.8X10-* 

Si 

39.9 

1.2 


*AU computed utuming a<Xp-iiux from observations, t).05/\]. 

^Si^silicata emiasivity, t/X» l/X emissivity, BB—blackbody (ci- constant). 


mtnimmng which is defined in the usual way. To 
allow for the instrumental effects, guiding noise, and 
other systematic errors, the uncertainties in the in^vid- 
ual points were constrained to be greater than or equal 
to a predetermined minimum error of 5% of the flux at 
that point The fitting algorithm was a grid search 
procedure adapted from Bevington (1969), which con- 
tinued until changed by less than 1% between itera- 
tions. It was found that the absorption feature in the 
SW spectrum was marginaily sif^icant statistically 
even though it is clear enough to the eye. 

The results of the fits are summarized in Table 4. 
The hot silicate/cold silicate fit is graphed in Figure 1; 
the other two types of fit produce curves which are 
nearly indistinguishable from that shown. It can be 
seen that the models give good fits to the spectra of the 
galactic center. This supports the plausibility of the two 
component model. 

V. 015CUSSI0N OF CONTINUUM MODELS 

All of the fits fail to produce as rapid a turnup at 16 
pm as observed. This can be ascribi^ to two causes: 
First, as mentioned above, a more realistic model would 
include a gradient of temperatures. Addition of a third 
component of higher temperature would tend to in- 
crease the 16 pm flux preferentially; but since the 
turnup is so sharp, it is unlikely that this is the main 
effect. The other possibility is that the true silicate 
emissivity is in fact more sharply peaked than assumed. 
In order to explain the discrepancy, the value of the 
emissivity at 16 pm would have to be reduced by a 
factor of ~0.8. 

Table 4 shows that the statistical significance of the 
fits is similar for all models, with the blackbody fits 
having slightly smaller x^- general, this means that 
we caimot say anything definite about the intrinsic 
spectrum of the hot dust or about its composition. A 


comparison of the NE and SW results shows that for a 
given type of fit the values of Ti, T^, and are nearly 
identical. Although the observed spectra are quite dif- 
ferent m both shape and normalizatiot , this difference 
appears to be due to the amount of hot dust in the 
source (i.e., f|), with the NE source containing 3 to 4 
times more hot dust than SW, Together with the simi- 
larity of t 2 , this suggests that the cold spectrum of SW 
might be better understood as due to the lack of a 
central heat source rather than an increase in mid- 
infrared dust extinction in front of an otherwise normal 
H a region, as suggested by Rieke, Telesco, and Harper 
(1978). The equilibrium temperature of dust grains at a 
projected distance of 30" from the NE source is esti- 
mated to be 30-60 K if the grain heating is due to 
radiation from the NE source. The SW source may be 
a local concentration of the “cold” dust near the galactic 
center. 

Note that the silicate models require about twice as 
much cold dust optical depth as the others. In the 
silicate models the hot dust cloud’s spectrum has a 
silicate emission feature, and more absorbing dust is 
required to hide it. Similarly, the silicate models require 
the highest optical depths of hot dust, while the fits 
the need somewhat less, and the blackbody fits need 
least of all, only about one-third of the silicate values. 
Finally, for all emissivity types, the cold dust is at 
temperatures near 50 K for the NE and SW sources 
and near 40 K for the 2'.7 beam. 

It is interestmg to extend this model to other wave- 
lengths to see if we can make it fail. To do so, we need 
a set of observations at comparable spatial resolution 
and covering a broad wavelength range. Luckily, the 
galactic center is a popular object, and Sgr A West 
(NE) has been examined at 30" resolution from 3.5 pm 
to 175 pm (Harvey, CampbeU, and Hoffman 1976; 
Sutton, Becklin, and Neugebauer 1974; Willner et al. 
1979). These data are plotted in Figure 4. Unfor- 
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Fio. 4.— Obfcrvatioiu of Sgr A Wot (NE) with a 30' beam, 
and modtii of the drut emiuion spectrum. Hatehtd arta, this 
paper. Cirelts, Harve^ *l aL (1970). Squans, Wilber *t aL (1979). 
Triatglta, Sutton *t aL (1974). Dot-ikuh lint, 8-13 pm spectrum 
of Woolf (1973), multlpUed by 0.1 for clarity, This spcctnim wu 
obtained with a 22' b^i and is shown here to deinonstrate the 
depth of the 9.7 pm sillcatfl fiuituie. Light solid iino, beat fit two 
component modri with both components consisting of silicate 
grains. Dottod lino, best fit two component model with hot dust 
having 1/X emissivity. Dashtd lino, three component model. See 
text for description. 


tuoately, due to experimental difficulties the 8-13 pm 
spectrum has not been measured with a 30* beam. The 
closest existing measurement is by Woolf and Gillett 
(Woolf 1973) with a 22* beam. It is plotted in Figure 4 
multiplied by 0.1 for clarity. The spectrum of Aitken, 
Jones, and Penman (1974) with a 13" beam and 15* 
N-S throw is similar in shape. However, Willner (1976) 
measured a number of compact sources at 10 pm with 
a 5* beam, including one spot O^beled ‘TRS 10") 
which was in fact mainly the 10 pm extended source. 
This spectrum shows a much shallower 9.7 pm silicate 
feature than the Woolf data. The lack of 30* data in 
the 10 pm window is particularly unfortunate because 
with it a direct comparison of the 9.7 and 18 pm 
silicate emissivities would be possible. 

The models are compared with the existing data in 
Figure 4. Here we have assumed the silicate emissivity 
used by Forrest, McCarthy, and Houck (1979) for 


A <16 pm. Only two of the two-component fits are 
shown;; the hot blackbody fit is similar to those plotted. 
All the models fit well in the 16-30 pm region, as they 
should, and all have similar failings os they ore ex- 
tended outside this range, 

Two serious problems with the model ore the lack of 
emission at shorter wavelengths and the lack of con- 
trast in the 9.7 pm silicate absorption feature. The 
former indicates that there are additional sources of 
emission at shorter wavelengths. The latter may indi- 
cate that the 9.7 pm opacity is larger relative to the 19 
pm opacity than has been assumed here. From the 
shape of the spectrum a temperature of about 400 K 
for the emitting material is indicated. A ratio of 
of about 2.5 rather than the 1.7 assumed 
is required to duplicate the depth of the 9.7 pm absorp- 
tion seen by Woolf. Such a ratio is consistent with 
values estimated from (S m] observations in obscured 
H II regions (McCarthy 1980). However, the com- 
plicated nature of the Sgr A source prevents a definite 
conclusion from being reached at this time. The data of 
Rieke, Telesco, and Harper (1978) and Becklin et aL 
(1978a) indicate the compact sources contribute an 
increasingly large fraction of the flux from Sgr A at 
wavelengths shorter than 20 pm. Therefore, to model 
Sgr A realistically at shorter wavelengths requires a 
knowledge of the intrinsic spectra of these sources. 

In the far-infrared, ail the two component models lie 
above the observations by factors of 1.2 to 3, but a 
simple and plausible explanation can be found. Harvey, 
Campbell, and Hoffmann (1976) fit their data with 60 
K grains having a emissivity and a 53 pm optical 
depth of 0.082. By using the law, t(19 pm)»0.58, 
or roughly one-t^d of the optical depth needed to 
explain the 18 pm feature in the silicate model (note, 
however, that r(53 pm) is in better agreement with the 
blackbody and 1/A fits]. This suggests that one should 
try three component models in which the new compo- 
nent acts only as an absorber. In order not to be 
observed in the far infrared, the dust must be either 
vety cold (T— 10 K) or extended (so chopping tech- 
niques would not detect it) or both. These properties 
are exactly what one would expect of interstellar grains 
lying far from hot stars, so we identify the third compo- 
nent with the interstellar extmetion in the galactic 
plane between the Sun and the galactic center. The 
second component is now interpreted as a local con- 
centration of cool dust (T~60 K) in the galactic center 
itself, and the ratio of column densities of local to 
line-of-sight dust is estimated at >—1:2, provided the 
hot dust grains are silicates. In Figure 4, a crude Et of 
this model is plotted, assuming T 2 ■•59 K, t 2 >0.6, and 
t 3 > 1.5. The parameters of the hot dust are the same as 
those of the two-component silicate model. 

On the basis of i^ared poiarimetry, Knacke and 
Capps (1977) arrived at a similar qualitative picture for 
the dust d^tribution. Observing that the 11.5 pm 



974 

polarization vector is roughly perpendicular to that at 
2.2 pm, they hypothesized that two dust components 
tvere involved. One is local to the galactic center and is 
responsible for the ll.S pm polarization. It is probably 
identical with our cool (T~60 K) local dust. On the 
other hand, Knacke and Capps suggested that the 2.2 
pm polarization is caused by another component dis- 
tributed along the line of sight. We identify this dust 
with our third, purely absorbing, component. 

The three component model is also supported by the 
similarity of the values of Tz for all three spectra. This 
fact suggests a screen rather than a shell of absorbmg 
material. Furthermore, the (S m] upper limits favor the 
models with hot silicate dust because of their higher 
dust optical depths. Finally, the far-infrared maps show 
a concentration of cold dust near the galactic center, 
and the existence of interstellar dust in the line of sight 
is well established. 

Since we have now derived estimates of t(19 pm), we 
are in a position to discuss the optical properties of the 
dust. The visual extinction to the late-type stars in the 
core has been estimated as —30 mag from the near- 
infrared colors (Becklin and Neugebauer 1968; Becklin 
er a/. 1978i). From the observed By flux at 2.17 pm, 
the extiactioa to the ionized gas is consistent with 
Ay -30 mag within large uncertainties (Neugebauer et 
at, 1978). However, several lines of evidence suggest 
that ..4 k h> not constant over the region (Lebofsky 1979; 
Rieke, Telesco, and Harper 1978), and we will consider 
Ay— 30 mag as a lower limit. From our models, we can 
set an upper limit of r(19 pm) ^ 2.1 for the central 30". 
Together, these give a lower bound of /4 k/t( 19 pm)^ 
14. Immediately this shows that the 18 pm feature is no 
stronger than the 9.7 pm one, since Ay/r(^.l pm) - 14 
toward VI Cyg No. 12 (Gillett et at. 1975). One might 
be tempted to derive a firmer number by assigning 
Ay -30 mag and r(19 pm)- 1.5 to the interstellar 
extinction, yielding ,4 k/t( 19 pm)— 20. This is risky, 
however, since the stars and gas appear to be em- 
bedded in the far infrared source. 

At this point, we can make an order of magnitude 
estimate of the dust mass column density in the line of 
sight This quantity, which we will call p, is related to r 
by p-t/kx, where is the optical ctsyss section per 
unit mass, as long as is constant along the line of 
sight The value of k\ is not well known, but at 9.7 pm 
it is thought to be —3X10^ cm^ g"' (Forrest, 
McCarthy, and Houck 1979). Since *^(9.7 pm)/(Cx(19 
pm)— 2, kx( 19 pm)— 1.5X 10^ cm^ g“‘. Therefore, the 
totd mass column density of dust toward the galactic 
center is .^1.4x 10 g cm~^, as long as both the hot 
and cold components are made of silicate grains. If the 
three component model is correct, then the mass col- 
umn density in the interstellar medium (ISM) is p(ISM) 
— 1.0X10 g cm~^, and the local component has 
p(local)— 0.4X 10~* g cm“*. At 10 kpc, our 30* beam 
covers an area of 1.6X10^^ cm^, so the dust mass 
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contained in this portion of Sgr A West is of order 3 
A/©. As shown by the 53 pm map of Harvey, Campbell, 
and Hoffman (1976), the cold local dust is in fact 
spread over a region roughly 90" X45" in extent, corre- 
sponding to a local dust mass of order 20 M©, 

VI, SUMMARY AND COKCLUStONlit 

In this paper we have presented threp 16-30 pm 
continuum spectra and two-color scans of the galactic 
center. Comparison with high spatial resolution maps 
shows that in our data the extended background emis- 
sion dominates over the compact sources. The map has 
a spatial resolution of 30" and shows that the source is 
sli^tly elongated along the galactic plane with a diam- 
eter of —40" or 2 pc. 

No [S in] 18,71 pm Us,e emission has been detected 
from Sgr A West. This is consistent with large dust 
extinction [r(18.71 pm) ^2,1] and normal S m/H n 
abundances. Indirect arguments suggest that S m/S ^ 
0.5 is reasonable, and that the S/H abundance is no 
more than a factor of 2 higher than in the solar 
neighborhood. A comparison with the By flux suggests 
a lower S/H abundance. Measurements of Ba in a 30" 
beam should allow more stringent limits on the S in 
abundance to be made. 

All three 16-30 pm continuum spectra show a broad 
depression centered at about 18 pm, which we attribute 
to absorption by cold silicate dust grains. In order to 
investigate the properties of the dust, we have con- 
structed two component models for the continuum and 
have fitted them to the data. In these modeb a hot 
optically thin cloud of dust lies behind a cold cloud of 
silicate dust grains with an optical depth of order unity. 
Three types of emissivities for the hot dust were consid- 
ered: a constant emissivity, a I/X enussivity, and a 
silicate emissivity. The derived dust parameters are 
luted in Table 4. All three types gave good fits, but the 
silicate type b the most plausible on physical grounds. 
The success of the models demonstrates that the as- 
sumed silicate enussivity b reasonable. Furthermore, 
the presence of the 18 pm feature supports the hy- 
pothesis that silicate minerals are indeed responsible 
for the 9.7 pm feature observed in the spectra of many 
objects. 

In extending this silicate model to X>30 pm, we 
find that a simple and natural modification improves 
the fit greatly. If one-third of the absorbing dust is local 
to the galactic center, and the rest b dbtributed along 
die line of sight (“interstellar”), then we can account 
for publbhed photometry to 175 pm. The temperature 
of the local dust is raised slightly, to —60 K, and the 
mterstellar component b too cold and/or too extended 
to be observed. 

Below X— 16 pm, the silicate model predicts a flux 
consistently below the observed values. Thu dis- 
crepancy requires the presence of additional emission 
sources at shorter wavelengths. 
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The silicate model also fails to duplicate the ob- 
served depth of the 9.7 /xm absorption feature. Tlie 
assumed value of 1,7 for the ratio t( 9.7 /im)/T(19 /xm) 
may be too low, A better value would lie between 2 
and 3, Finally, a lower limit of 14 may be set on the 
ratio of visual to 19 /xm extinction, Wf'/r(19 /xm). 

In the line of sight to the galactic center, the total 
dust mass column density is —l.dX 10"’ g cm"’. This 
value assumes that silicate grains make up both compo- 
nents; other models will lower this estimate by a factor 
of 2 or 3. In the three component model, the interstellar 


dust mass column density is —I. OX 10”’ g cm"’ along 
the 10 kpc line of sight. The lofal dust mass I ;cal to the 
galactic center (within about 2 pc) Is —20 Mq, 
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